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The x-ray interference is calculated for layer lattices in which the phase shifts between con- 
secutive layers and the scattering powers of individual layers do not follow a strictly periodic 
arrangement. In the second section the scattering power of all layers is assumed to be the same 
but the phase shifts can take on different values. In the third section neither the scattering 
powers nor the phase shifts have fixed values but a simplifying assumption is made about the 
; phase shifts according to which distances between neighboring layers can be represented as 
sums of two distances characteristic of the individual layers. In both these sections a random 
sequence of the layers is assumed. In the fourth section the problem of arbitrary scattering 
powers and phase shifts is treated, and furthermore a statistical correlation between neighboring 
the layers is introduced. In the following section the general theory is applied to a specific partially 
_— ordered stacking of layers encountered in micas and other similar minerals. The last section 
on treats irregularities in close packed structures of spheres and irregular sequences of layers in 
—_ graphite. 
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smoothly and monotonously with angle the 
former shows sharp discrete maxima due to the 
infinite extension of the layers in two directions. 
Layer form factors have been calculated by adding 
the scattering from the atoms within the layers.' 

The amplitude of the radiation scattered by 
the crystal may be obtained by adding the am- 
plitudes scattered by the individual layers. This 
reduces the calculation to a single summation 
over layer form factors and thus leads to a one- 
dimensional problem. The mathematical methods 
used can equally well be applied to the calcula- 
tion of reflections from a single layer consisting 
of parallel but otherwise not strictly regular 
sequence of atomic chains. Again the scattering 
from a chain with irregular sequence of atoms 
could be calculated in the same manner. Prac- 
tical applications, however, are limited to layer 
lattices. 

Only problems will be treated in which there is 
some phase relation between scattering from dif- 
ferent layers. If no such phase relation exists the 
total intensity is obtained by adding the inten- 
sities from the individual layers.” 

In the second section it is assumed that all the 
layers have the same form factor. The phase 
shifts of radiation scattered by neighboring 
layers are supposed to have a limited number of 
values. Each phase shift occurs with a certain 
probability but there is no correlation between 
consecutive phase shifts. A practical example for 
which the above assumptions are approximately 
valid is the x-ray scattering from the hydrated 
clay mineral montmorillonite.’ In this mineral 
identical layers containing heavy atoms are 
separated by varying numbers of water layers of 
considerably smaller scattering power. If we 
assume that the number of adjacent water layers 
is limited and if we neglect the scattering by 
water layers the influence of the water will be to 
produce a limited number of phase shifts between 
the scattering waves of neighboring heavy layers. 
The calculation yields intensity fluctuations 
rather than sharp interference maxima _ cor- 
responding to orders of reflection. 

In the third section crystals are considered 
which consist of layers with different form 


1M. von Laue, Zeits. f. Krist. 82, 127 (1932). 

2B. E. Warren, Phys. Rev. 59, 693 (1941). 

3S. B. Hendricks, R. A. Nelson, and L. T. Alexander, 
J. Am. Chem. Soc. 2, 1457 (1940). 
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factors. It is assumed that each type of layer 
has a definite probability of occurrence, but the 
sequence of layers is random. Distances between 
neighboring layers and correspondingly their 
phase shifts depend upon the kind of layer. In 
order to obtain an explicit expression for the x-ray 
interference it is further assumed that the 
distance between two neighboring layers is a 
sum of two lengths characteristic of the layers. 
This assumption is satisfied if crystal distances 
can be obtained from ionic radii. A practical 
example is the precipitation of copper in a 
duralumin alloy in which the copper segregates 
on (100) planes of the crystal. Another example 
is afforded by the basic salts of the divalent 
metals, for instance n CdCl.-m Cd(OH)s. Here 
CdCl. and Cd(OH). layers may follow each 
other in a random sequence. 

In the fourth section a more general problem 
is treated than in the third. It is no longer 
assumed that distances between neighboring 
layers are sums of two characteristic lengths. 
Furthermore, a correlation between occurrence of 
layers is introduced such that the probability of 
a layer to be present in a certain position depends 
upon the neighboring layers as well as on the 
abundance of the layer in question. This cor- 
relation affords a transition between the case 
of a random sequence and a completely ordered 
crystal. This problem is not solved explicitly but 
is only given in terms of the solution of a secular 
equation of the mth order where n is the number 
of different kinds of layers. The practical 
examples are the same as given in the preceding 
sections only that a more exact and more labo- 
rious solution is offered. 

In the final two sections the secular deter- 
minants are solved and intensities of interference 


explicitly given for two cases to which the treat- 


ment of Section 4 applies. The first of these is 
the stacking of layers in micas, kaolinite, etc. All 
layers have the same form factor but they may 
occupy any of three positions differing by shifts 
parallel to the plane of the layer. This is possible 
since relative positions of neighboring atoms in 
the top of one layer and the bottom of another 
are not affected by such shifts. Probabilities are 
assumed for each possible relative shift of neigh- 
boring layers. The second example treats irregu- 
larities in close packed structures of spheres and 
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irregular sequences of layers in graphite. In both 
these cases regular packing obtains for neigh- 
boring layers and irregularities enter only in 
relationships of second neighbors. Nevertheless 
the problem can be reduced to the one treated 
in Section 4. Neither of these cases has yet been 
observed’ but the occurrence of the graphite case 
is to be expected since there is evidence for per- 
fectly random shifts parallel to the planes of the 
layers in some crystals. 

Influences of finite size of crystallites, thermal 
motion, divergence of x-ray beams, and the 
dynamical interference theory are not considered 
in this paper. 
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Two publications by Landau and by Lifschitz‘ 
give explicit intensity formulae for several types 
of arrangements of layers assuming no correlation 
between the layers. The following Section 2 of 
the present paper gives a somewhat more general 
treatment of the same problem. The influence of 
the correlation between layers has been con- 
sidered later by Lifschitz’ for the special case 
that the crystal consists of two kinds of layers. 
The calculations given in the first part of Section 
4 of the present paper deal with the same 
problem. These calculations have been included 
merely as an introduction to the treatment of 
crystals containing several kinds of layers. 


SECTION 2 


It is assumed that all layers have the same form factor. For the sake of simplicity the case where 
only two phase shifts, g1, and g2, may occur between neighboring layers will be considered first. The 
frequency of occurrence of the two phase shifts is assumed to be equal and the sequence of the phase 
shifts is random. The alternation of phase shifts between two values may be due either to displace- 
ments of the layers parallel or perpendicular to their planes. 

The amplitude and phase of the scattered radiation is represented by a vector in the complex plane 


\V| =| V| exp ig, 


where V is the amplitude of the wave and ¢ stands for its phase. If V, represents the scattering 
from the kth plane the scattering from a crystal consisting of n planes is 


De V:.. 
1 


The absolute square of this sum gives the intensity. The intensity per layer is 


1 


n 


T=-DVkRDV:- Vi, 
1 1 


where V,* is the conjugate complex of V;. If m approaches infinity one may replace the double sum- 
mation and subsequent division with n by a single summation in which the product of one vector 
with every other vector of the lattice appears. Thus 


+2 +x“ 
I= (3) DLA(ViVo*+V.AV0) = DR | V;. | Vo| cos (ge— Yo). (1) 


Here (g,— go) may be written as 


(Gx — Go) = (Gk — Gr—1) #(Gr—-1 — Gr—2) H+ (G1 — Go). 


Each of these differences may have the values ¢ or ¢® with equal probability, where g“? and ¢®) 
are the two possible phase shifts between neighboring layers. 


m-___——____. 


*L. Landau, Phys. Zeits. Soviet Union 12, 579 (1937); I. M. Lifschitz ibid. 12, 623 (1937). 


*I. M. Lifschitz, J. Exper. Theor. Physics U. S. S. R., p. 500 (1939). 
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The average intensity per layer I, is obtained by averaging Eq. (3). Since the absolute values of 
the vectors are all equal the actual averaging need be carried out only over the cosines. 


+2 +0 
In = DRE! Vi! | Vol cos (ge— Go) w= V? SRLcos (yx— G0) Iw. (2) 
bo bao 
Here V? is the absolute square of the scattered vector for any layer, and the average over the sum 
has been replaced by the sum over the average values. Since 


[cos (go— go) Jw=1, [cos (¢:— 0) Jw= [cos (¢-«— 0) Jw, (3) 


the calculation need be carried out only for positive k values. Without loss of generality the phase 
of the scattered wave can be normalized in such a way that go=0. If 3(¢@+¢°) =@ it follows that 


[cos gx Jw=[cos ké cos (y—kS) —sin kG sin (g,—RS) |wv 
=cos kg[_cos (¢:—kd) Jw —sin kG[sin (¢,—kS) Jn. 


The arguments (g;—k%) may assume values differing in sign only and such a pair occur with equal 
probabilities so that [sin (¢.—k¢) ]4=0. Thus 


[cos gx: lw=cos RGLcos (g,—RG) Iw. (4) 
The average values of the cosines can now be found by complete induction. 
Loos (¢i41—(R+1) 4) Iw=[Leos (yg —kG) cos [(ge41— ge) — Sw 
—[sin (¢-—R@) sin [(grs1—¢r) —S] Jw. (5) 


Since there is no correlation between subsequent phase shifts the average over the factors of the 
right side may be carried out independently. The average over the sin terms vanish. According to 
the definition of @ 

(gi41— Gr) — = £3( 9 — yp) 
and by complete induction 
Leos (¢ — kG) Jw= (cos 3(¢ — o))*. 


From (2), (3) and (4) it follows that 
Iy= V?{1+2 rk cos kg(cos $(g" — v))*] 
= V>L1+ZE(cos Ho" — ¢))* exp ikB+ Lh (00s He — o))* exp ike) 
Evaluating the geometrical sums and writing the results with a common denominator J,, is obtained as 


(6) 





I r| 1—[cos 3(g) — g®) 
i= | 


1—2 cos $ cos 3(g — ¢) +[cos 3(g? — o) ? 


This formula gives in general a scattered intensity in every direction for which V is different from 
zero. The influence of the layer form factor, of course, is contained in V. If both g“ and g® approach 
multiples of 27 at the same rate as their difference approaches 27 then Eq. (6) becomes infinite. In 
fact while the numerator approaches zero linearly the denominator approaches zero quadratically. 
This corresponds to a sharp interference maximum. 

In order to illustrate intensity fluctuations given by Eq. (6) scattering of a layer lattice in which 
the distances of neighboring layers have the values 10A and 13A with equal probabilities will be 
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Fic. 1. Scattering from a random mixture of two layers having equal scatter- 
ing power but different phase shifts corresponding to interlayer spacings of 
10A and 13A. Angle of incidence equals one-half of scattering angle, 6. Results 
according to Eq. (6). 


considered. This lattice corresponds to the partially hydrated montmorillonite where 10A obtains 
for the anhydrous material and 13A obtains if there is one water layer between two silicate layers. 
Figure 1 gives the intensity of reflections from the parallel set of planes as a function of (sin/A)6/2 
where 6/2 is the angle of incidence. The full perpendicular lines correspond to consecutive orders of 
reflections for a lattice with layers having 13A as their separation, while the dotted perpendicular 
lines refer to a lattice with a 10A constant. It is interesting to observe that the curve representing 
Eq. (6) shows marked maxima where sharp interference of these two crystals occurs near the same 
angle while lower and broader maxima are found for isolated interference maxima of the regular 
lattices. Of course no definite orders can be assigned to the intensity fluctuations in the curve. For 
(sin/A)6/2=0.5A~ the 13th order for the lattice with 13A separation and the 10th order for the 10A 
lattice coincide and at this point Eq. (6) becomes infinite corresponding to a sharp interference 
maximum. The actual substance seems to follow the behavior shown in Fig. 1. 

In generalizing the above results it is assumed that r different phase shifts g“, g®, ---, g can 
occur betwee neighboring layers with frequencies f“, f°’, ---, f°. Since one of these phase shifts 
must be presett between any two neighboring layers the sum of frequencies is 

\ 







ds fP=1. 
1 


It will again be necessary to calculate Eq. (2) which amounts to calculating [cos (g— go) Jw. 
Eq. (1) remains valid but each term on the right-hand side might have any of the values ¢“) (s=1, 
2, --- vr). Essentially the same reasoning can be used as before. An average phase shift 2 is defined 
by the implicit equation 


’ 


Ds f® sin (g® — g) =0. (7) 
1 


It is necessary to show that [sin (¢g,—k¢) ]4=0, a relation which in the previous case followed 
immediately from equal frequency of occurrence of positive and negative values of the argument 
for the sine. One finds 


[sin (G41 —(R+1)8) Jw=[sin [(¢r41— or) — S] cos (gr —kG) +008 [( G41 — Ge) — S] sin (yr — RS) Jn 
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The averaging can again be carried out independently over the factors. Since (¢+1— ¢.) must have 
one of the values ¢“? with the frequency f the average of the first factor in the first term on the 
right-hand side vanishes according to Eq. (7). Thus 





Csin (¢n41—(k+1)$) Jw=Leos ((¢e41— gx) — 2] sin (gr —kO] Jw 


and it is apparent that [sin (g4:—(k+1)6=0 Jw if [sin (g,—k@) Jw=O0. Since Eq. (7) means that 
[sin (g,—k¢)=0]4 for R=1, one obtains by complete induction that [sin (¢,—k%)=0], for any 
value of k. 

Now Eg. (4) and a relation 


[cos (PK4-1 _ (R+1)d) av = [cos ( gx —ko) |wLcos C(eepi— $x) sity % Jn, 





which is analogous to Eq. (5) follow in the same way as in the previous proof. 
Introducing the abbreviation 


Loos ((¢rs1— gx) — 4) w= Ls f cos (g — g) =C, (8) 
1 


“it can be shown by complete induction that 
[cos (¢:—kR%) Jw=C* and [cos yy Jw=C* cos kG. 


Equation (2) now assumes the form 
V2L1—C?] 
7 1—2C cos o+C? 





AV 


This expression is the same as Eq. (6) only the quantity cos (3)(¢“’ — ¢®) has been replaced by C 
and the quantity ¢ is defined by Eq. (7). 


SECTION 3 


A mixture of 7 kinds of layers differing in their form factors is now considered. The sequence of 
layers is random and, at first, it is assumed that the phase shift is the same between any two con- 
secutive layers. 

Considering a finite crystal of m layers the radiation scattered by the kth layer can be written as 
V“ exp ikg; here V“ is the scattering vector descriptive for the rth kind of layer which is supposed 
to occupy the th position. The quantity V“) depends of course on the directions and polarizations 
of the incident and scattered x-rays beams but it does not depend on the phase shift between layers. 
This latter dependence is contained in the factor exp ikg. It must be borne in mind that the V“? are 
in general complex, for instance translation within the plane of a layer multiplies the V“ by complex 
factors with the absolute value unity. The average intensity per layer can be written as 


1 n n 
In=~— Dk LI VV* exp i(k—De, 
nN 1 1 


where V“)* is the scattering of the ¢th kind of layer which happens to occupy the /th position in the 
sequence. The averaging is to be performed over (s) and (¢) each of which can independently take { 
on values from 1 to r with frequencies f™, ---, f‘ which again satisfy : 


ds fM=1. 
1 
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Thus the average intensity is 


1 n r 1 n r 
In=— X/RI LstlVOVO*f fO exp i(k-D gt+— Xk Ds FO! VO !27, (10) 
n 1 1 n 1 1 


where the prime on the first summation indicates that k=/ must be excluded. The occurrence of two 

sums is due to the fact that for k#/, f* and f are independent probabilities while for k=] both 

s and ¢ describe the same layer. Since in the second sum on the right-hand side k does not occur under 

the summation sign, the summation over k giving simply a factor n which cancels the factor 1/n. 
The scattering vectors V“) can be consider: d as sums of two vectors 


Vv =U4+W™, (11) 
where U is an average scattering vector defined by 
U=d's fVv™. 
1 


The W“ are then defined by Eq. (10). Multiplying both sides of Eq. (10) with f and adding 


r 


> (s) SOW =() (12) 
1 
is obtained. 


Substituting (10) into (9) shows that the average intensity can be written as a sum of three terms 


Ty _ Laut Lut | 


where 
1 n r r 
lun =— DRI Sstf © fOUU* exp i(k-D ot Es f! U2, 
1 1 1 


n 


ji— 


Law =- dhl Vstff@TTW*4+wour) exp i(k—l)ot+ ds fO(TW*+Wou*), 
1 1 1 


= 


Dew =— DRL Dstf © f OWOW* exp i(k—I e+ Csf| WO |2. 
n 1 1 1 


In J, the summation over s and ¢ gives according to Eq. (9) 


1 n» n 
Tin=— Lk YI) U\2 exp i(k—De (k=I not excluded), 
1 


nN 1 


which is the same result as would be obtained for a completely regular layer lattice where all layers 
have the form factor U. The summation over s and ¢ shows that the expression J,,,. vanishes according 
to Eq. (12). The same holds for the first sum in J,,.. while the term 


$s f@| W® 2 
1 


gives the average scattering of a single layer having scattering vectors W“ with probabilities f. 
Thus the two non-vanishing terms in J, correspond to the scattering from a completely regular 
lattice and to scattering uninfluenced by phase relations in which the angular dependence of the 
scattering is determined by the layer form factor alone. In other words the crystal behaves as though 
it were completely regular and had a form factor equal to the average of form factors and as though 
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there were present in addition layers scattering independently and having form factors equal to 
the deviation of each real form factor from the average. If only two kinds of layers are present it 
can easily be seen that in the latter kind of scattering the difference of form factors enters. 

The results here obtained can readily be generalized to a random mixture of scattering centers 
in any kind of lattice. In fact phase relations do not enter explicitly in the argument and as long as 
phase relations are the same for all r kinds of scattering centers the above result remains valid. The 
results for the general case have actually been found by von Laue in 1918.° Thus the scattering from 
a completely random solid solution gives interference corresponding to the completely ordered 
crystal with an average form factor on which a continuum is superimposed with an intensity dis- 
tribution that depends on the deviation of the individual atom form factors from the average. This 
continuum will have an appreciable value at zero angle. 

An explicit formula for the scattering can be obtained even if it is assumed that the phase shift 
between two neighboring layers of the kinds s and ¢ depends on the nature of these layers. But, a 
simple solution can be obtained only if this phase shift ¢“" is the sum of two phase shifts g® and 
¢‘® characteristic of the separate layers; g“ and ¢“ will of course depend on angle of scattering as 
well as the wave-length of x-rays. The average intensity then becomes 


1 n n 
n=|- Dk YI V.OV* exp ieu-| ; (13) 
1 


n 1 Av 


where the subscripts k and / emphasize that the vector V;.“? is in the kth and the vector V,“ in the 
ith position. 

The angle g.—,) is the phase shift between the &th and /th layer. Since gi—1, depends not only on 
the nature of the kth and /th layer but also on the nature of the intervening layers the averaging 
cannot be replaced now by summations over the two layers involved but complete induction must 
be used as in the previous section. 

Equation (13) can be transformed in the same way as the expression for the intensity was trans- 
formed in the beginning of Section 2. Writing | V*| =| V*| exp 7a“ one obtains for 


+x 
In= Dk V, 1 Vo cos (gu—0y ta —a) |py. (14) 


The summation now extends from — x to +< the crystal having been replaced as in the second 
section by an infinite sequence of layers. ¢%—9) is the phase shift between the kth and zeroth layers. 

The average over the kth term in Eq. (14) will first be considered. Assume k>0 and start with the 
averaging over the (k—1)th layer by substituting all the possible kinds of layers into that position 
but leaving the other layers unchanged. Assuming that the (k—1)th layer is of the pth kind and the 
(k—2)nd layer is of the gth kind one may write 


PY (k—0) = P(h-2)-0 FE +29 + . 


For the average value over the (k—1)th layer is obtained 
yp f® | Vi | | Vo! cos Cecu—2y)- HEM +29 + ge +a —al ] 
1 
‘ 
ia <p f@ V;. | | Vo" | cos [ ecu—2)-0 HE + 9 +a” —al +3] cos (2y sie. 2) 
1 


° 
_— Pp f™ | Vi || Vo“ | sin [ecu—2)-0 HE + 9° +a —a% +o] sin (2y a Q). (15) 
1 





6 M. von Lane, Ann. d. Physik 56, 497 (1918). 
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If in this expression ¢ is defined by the implicit relation 
r 
Dp f™ sin (2¢™ — g) =0, (16) 
1 
the second term on the right-hand side of Eq. (15) vanishes. Introducing the abbreviation 
q 
. 
> p f™ cos (26 — 3) =C (17) 
1 


and performing the average over the remaining layers which have been kept fixed so far there results 
[| Ve] | Vo| cos (ga 0 ta —a) Jay 
= [C| V,.@ | | VV," | cos (O¢(k—2)—0) a got yg) +a) —ai) + 2) Inv (18) 


Since g“ + ¢ would be the phase shift between a layer of the gth and sth kind if these layers were 
neighbors the factor by which C is multiplied in Eq. (18) differs from the average taken over a pair 
of layers in the positions 0 and k—1 only by the phase ¢ appearing in the cosine. One can, therefore, 
write 


[| Vi. | | Vo" | cos (¢¢:—0) ta —a) Ja = [C| V (1) || Vo cos (¢e—1)-0) fa” —a+ 6) |p. 
By complete induction it finally follows that 


[| V;,.@ | | Vo | COs (ou ota®™—a™) |y=LCe-Y V1 Vo cos Yul ota —a +(kR—-1)d) |av 
r r 
=Ps Vt CO-VFOFO! Vy! | Ve] cos (9+ e%+a™—al+(k—1)8). 
1 1 


For negative k values the same expression is obtained except that |k| rather than k appears in the 
formula. The contribution of the term with k=0 to Eq. (14) is 
r 


Ls fOlVO|2, 


1 


Thus is obtained 


7 


In= Ls fO| VO |2$2 Dk VstC#-V FO FO | VO| | VO| cos (G+ ge +a —a +(k—1)B). (19) 
1 1 1 


The cosine in the second term in Eq. (19) can be expressed in exponentials and the summation is 
then reduced to a geometrical series giving 


r r 
L.= Ds fO|VO|242 Yst fOfe| VO| | VO| 
1 1 


cos (¢@ + g +a —a) —C cos (y+ op — +a —a 


icsiaesasegiaabie -?. (20) 
1—2c cos 6+C? 





It is of interest to compare this formula with the one obtained in the first part of this section. In 
this case the ¢) become equal to each other and to ¢/2. From Eq. (17) it follows that C=1. After 
breaking down V“ into U and W®) Eq. (20) can be transformed into two terms the one containing 
the W’s and giving no interferences and the other containing U and giving sharp interference 
maxima. 
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SECTION 4 


The case where the sequence of layers is no 
longer random, but where correlation between 
neighboring layers is introduced will now be 
considered. To illustrate the procedure to be fol- 
lowed in this case it will first be assumed that 
two kinds of layers occur with equal frequency 
and that the phase shift between neighboring 
layers is always the same. The layer form 
factors are real and have the values V™ and V®. 

To introduce correlation between layers two 
neighboring layers are considered and it is 
assumed that the probabilities that a layer of the 
first kind is followed by a layer of the first kind 
is A,“), for a layer of the first kind by second 
kind A,“*), for a layer of second by first kind 
A,®), and for second kind by second kind A,?. 
Since the first layer of the pair has a probability 
of 3 to be of the first kind and 3} to be of the 
second kind it follows that 


A,OD+4A,02 =F, A,@D+4A4 22) =F, 


Since the second layer of the pair has also the 
probability 3 to be of the first or second kind one 
finds 


A,{QD+A4 2) =i and A {244 1022) = 3, 
From these relations it follows that 
A\O®=A,@0 and A,OD=A,22, 
The following abbreviation is now introduced, 
AOD =A,2) =144, 
A,O%) =A,20 =1— A, 


The probability P“ =2A“ that the second 
layer of a pair is of the kind one if the first is of 
the kind one is now introduced. Similarly 


POX =2A02%, PeAD=2A42), P22) =724 (22), 


(21) 


Then the correlations between second neighbors 
can be calculated by multiplying the P’s. In fact, 
the probability for a third layer to be of the first 
kind if the first was of the first kind is 


POD = P,Av. Pad 4 P.O» Pew 


and the probability A»“" to find two layers of the 
first kind as second neighbors is given by 


AOD =1P,0) = 274,00. 4,004 4,02. 4,207. 
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Similarly the probabilities of other combinations 
as second neighbors are 


A,°2) =  P{22) - 2[A 1°22) -A,CX+4A4 20 -A 12) 7], 
Ag!) =3 P02) = 20 A OD. 4 0944 02). 4 227, 
Ag2@) _ 3 P, (20) = 2[A 122) -A ;2D+A4 yb -A yD 7, 
From Eq. (21) it follows that 
Ao) = A o(22) =1[1+(4A)?*], 
Ao” = Aol?) = 4 = (4A)?)]. 
In a similar manner the probabilities for any two 
kinds of layers occurring as kth neighbors, A;,“), 
A, A;,“@, A,@ can be introduced. These 
probabilities can be reduced to the probabilities 


of finding 2 kinds of neighbors as (k — 1)th layers 
by the relation 


AOD =2(AO Gy) +A OD +A O® G1) +A @P) 


and three similar equations. One can easily 
verify that the equations for arbitrarily removed 
neighbors are satisfied by writing, 


A,O) = A,@2 =1(1+(4A)*, 


(22) 
A,@Y = A, =1(1—(4A)*). 


For k=0 we have 


Ay =Ay® =4; Ag?) =Ay@ =0 


which may be interpreted as meaning that a 
layer cannot be of the kind 1 and 2 at the same 
time and has the probability } of being either 
kind 1 or 2. Thus Eq. (22) holds even for k=0. 

The average scattering per layer which accord- 
ing to Eq. (1) is 


+00 
In = UR 2(ViV0*+Vi*Vo) 


is now calculated. Since the phase shift between 
neighboring layers has been assumed to have a 


. fixed value, ¢, V; is then given by V™ exp zk¢ or 


V®) exp ikg according to whether the kth layer 
is of the first or second kind and Vp is simply 
V™ or V®. Then the average intensity is 


+00 Tr 
Inw= Dk Sst Avo? VV cos ke. 


—2 1 


The averaging has been replaced by multiplica- 











= {{_ —_— —- 


ons 


|, 
I, 
)'}. 


two 
k (11) ; 
hese 
ities 
yers 


)) 


asily 
oved 


(22) 


lat a 
same 
ither 
k=0. 
cord- 


tween 
ave a 
ky or 

layer 
imply 


iplica- 








PARTIALLY ORDERED LATTICES 157 


tion with the probability factors A,“ and sum- 
mation over s and ¢ each of which can assume the 
values 1 or 2. Using Eq. (22) it follows that 


in VOLVO 2 
n= EH ( —) cos ke| 


‘aii VM — pe: 
+EA ( ; ) (4A )* cos ke 











The first term on the right-hand side corresponds 
to the scattering of a perfectly regular layer 
lattice in which the layer form factor is 
(V+ V®)/2. The second term gives an intensity 


yo— V® 2 1—(4A)? 
| = eg eee = . 
( 2 1—2(4A) cos ¢+(4A)? 





(23) 


If A=0 we have complete disorder and the term 
reduces to [(V“ — V®) /2 ? which is the result 
found in the first part of the third section. A is 
positive corresponding to a greater probability 
of like layers to follow each other. Maxima of 
intensity occur if g=n"X2z, i.e., in the same 
positions where the regular lattice has its sharp 
maxima. If A attains its maximum value } the 
crystal breaks up into regions containing ex- 
clusively one or the other kind of layer. In this 
case the second factor in Eq. (23) becomes zero 
unless ¢ is »X2m when both numerator and 
denominator vanish. For such values of ¢ the 
denominator approaches zero quadratically if 
(44) approaches unity while the numerator 
vanishes linearly so that for this ordered case 
Eq. (23) gives sharp interference maxima the 
intensities of which are proportional to 


[(V™—Vve)/2¥. 


These maxima have the same positions as in an 
ideal crystal and the sums of the intensities 


C((VO+4V)/2P4+[(V —V@/29 


give the average intensity of two crystals with 
the layer form factors V™ and V®). 

For negative values of A Eq. (23) has its 
maxima if g is an odd multiple of x. This gives 
half orders lying between the maxima of the 
regular lattice corresponding to the fact that for 
negative A values there is a tendency for opposite 
layers to follow each other. These maxima be- 


come sharp for A = —} when the random lattice 
degenerates into a regular superlattice with the 
layers of the first and second kind alternating. 

The most general case is now considered. There 
shall be r kinds of layers possessing complex scat- 
tering vectors V, ---, V“ and occurring with 
frequencies f“, ---, f°. The phase shift between 
two consecutive layers of the s and ¢th kinds 
shall be ¢“”. It will frequently be true that 
gy") = 9‘) although this holds necessarily only 
if the two sides of each layer are equivalent. 
There also shall be correlations between the 
probabilities of occurrences of neighboring layers. 
But the only correlations of probabilities between 
second and farther neighbors shall be those that 
follow from correlations of immediate neighbors. 
An explicit formula for the intensities of scattered 
x-rays cannot be given. The final formula will 
contain expressions that depend on solutions of 
a secular determinant of the rth order. 

Denote by A,“ the probability that a certain 
layer shall be of the kind s and the following 
layer, in a direction considered as positive, of the 
kind ¢t. The sum 


. 
>s A,“ 
1 


gives the probability that the second layer is of 
the kind ¢ and therefore must be equal to f™. 
One finds in a like manner that 


Ft Aye) =f, 
1 


Similarly A.“ gives the probability for a layer 
of the kind s to be followed by a second nearest 
neighbor of the type ¢. In general A,“*” gives the 
probability that a certain layer is s and that its 
kth neighbor in the positive direction is ¢. It 
again will be frequently true that 


A,@9 =A, (24) 


Also the matrix A; with the elements A,“ is 
introduced. The matrix is real and if Eq. (24) is 
satisfied it is symmetrical. 

Just as in the special case at the beginning of 
the section it is useful to introduce numbers 
P,“ which give the probabilities that the kth 
neighbor of a given layer shall be of the ¢th kind 
if the layer itself is of the sth kind. These numbers 
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are connected with the A;“° by the relation of the s kind. And the sum of these probabilities 
A,@) = f@ Po, over all possible values of q will give the total 
probability of s being followed by tas its (k+1)th 


Matrices with elements P;,“® will be repre- neighbor. 
sented by P,. It is then easy to see that the From the above statement about multiplica- 
matrix product of P; and P; gives Py.41) because tion of P matrices it follows by complete induc- 
the probability P¢.41)°° may be written as tion that 
P,.=(P)* 
P 
st) —. (84 t) ry. . . . . 
P esp = og POOP. The diagonal matrix F with the diagonal cle- 
1 


ments f°) =f is also introduced. Then 

The sum on the right-hand side gives the prob- . 
2 . . ing A, = FP, = F(P))'. 
abilities of the arrangements of g and ¢ layers in 


the k and k+1 positions if the original layer was The average intensity per layer can be written 





In = rk VstL A. 23 (VoOV.OF+-V_*V ©) Jay 
—2 1 





= 2s f®| Vo 24k Ds tl Ap ©VoV.* Jat vk dst [A pO OV M*V, Jy, (25) 
1 1 1 1 1 


+o 
The latter expression arises from breaking the sum 0k into a term for k=0 and two others 
is) +o —@ 
for >-k and >¢k. The last two terms in Eq. (25) are conjugate complex and the last can therefore 
1 -1 
a) 


be written }>k over conjugate complex terms. The scattering vector Vo“) may be set equal to the 
1 
layer form factor V“? while V,“° is the product of the form factor V“ and the factor exp i¢,. Here 
exp 1g," is the phase shift between the zeroth and kth layers which are of the sth and ¢th kinds. 
This phase shift depends of course not only on s and ¢ but also on the kind of the intervening layers. 
The averaging in Eq. (25) must be carried out over these intervening layers since the averaging over 
the zeroth and kth layers has already been included explicitly in the summation. 
Thus Eq. (25) may be written as 


In= ds f | Vo | 2+ vk DL stVV*A,©9L exp (—t¢.@?) wt DR Dost VORVOA.OLexp ig? Ja, 
1 l 1 1 1 


where the averaging need be carried out only over the exponential factor. This averaging can be 
carried out by introducing further matrices Q; with elements P;.“ exp (—ig,“”). Then the relation 
holds 


Q41) = QQ. (206) 
This can be seen in the following way 
. 
(st) = p*”  ¢ 2 0 = P,P exp (—i( gy + (@)) (27) 
Qos =P usr pLlexp (—t¢ e+) Jw= Lg Pi pe EXP (— UA Pk g)). - 
: 


Indeed Q“" ¢.41) is the sum of the phase shifts between layers of s and g kinds occupying the zeroth 
and kth positions plus the phase shift between two neighboring layers of the gth and ‘th kinds which 
have been denoted by 9%”. 

Since averaging over g is explicitly included in the summation the averaging process need be 











al 
un 


ers 
ore 
the 
ere 
ds. 


ers. 
ver 


Jas 


1 be 
tion 


(26) 


roth 


‘hich 


d be 











PARTIALLY ORDERED LATTICES 
carried out only over the factor exp (—ig,.“”) so that Eq. (27) can be written 


6) 2 a r 
Qas+y = Lg Pr Lexp (—t¢.%”) JP i? exp (—719%) = Vig Q.PQ 
1 1 
which proves that Eq. (26) holds. 
It follows from Eq. (26) that 
O. = (Q)". (28) 
With this notation 74 may be written 


r 


y= 72 f | V&) | 24 dk dst VOrFVO<IOO,.COLSR dest VOV orp, @0* 
1 1 1 1 


1 
=>s f™| V@| 245k spur VFQ,.+ > k spur V*FQ,*. (29) 
1 1 1 


Here the symbol spur means the diagonal sum of the matrix product. Matrix Q,* is a matrix whose 
elements are conjugate complexes of the elements of matrix Q;. The matrices lV and !’* have the /, 
s elements V°O*V™ and V‘0V™*, respectively. 
x 
In order to carry out }°k it is convenient to transform Q, to the diagonal form. One can show 
1 
that this is always possible if one of the two following conditions is satisfied ; all proper values of Q, 
are different; the matrix A is symmetrical and ¢“” = ¢“*). (In the latter case the transformation of 
Q to the diagonal form can be reduced to the transformation of an Hermitian matrix to its diagonal 
form which can always be carried out.) If none of these conditions holds then it may not be possible 
to diagonalize Q). and a formalism is needed that is somewhat more complicated than the one to be 
presented here. 

One may diagonalize the Q matrices by introducing a matrix O for which OQ,O~ is diagonal. The 
diagonal elements of this matrix, which are denoted by Q“?, are the solutions of the secular equation 
belonging to the matrix Q;. The matrix elements of O may be obtained from a system of homogeneous 
linear equations whose determinant is the secular equation. The same transformation matrix O 
that diagonalizes Q; also diagonalizes Q;. This is true because of Eq. (28). 

The diagonal elements of OQ;.0~ are the quantities (Q)*. The spur occurring in Eq. (29) may now 
be written as 


spur V FQ; =spur (OV FO™)(OQ,0~'). 


The second factor is in the diagonal form and therefore only the diagonal elements of the first factor 
enter into the spur. If the first factor is 
R=OVFO"! 


whose diagonal elements are R“*) one obtains for In, 
r 4) r i 4] r 

n= Zs fO| VOP?ALR Ys ROMO) TR Ds RQ) 
1 1 1 1 1 


r r ROOQ® - ROO*OC * 
= dos f®! V@|24+ 5s ———-+. }'s —___- 
7 I i ant Qe 1 = Qe 
The first term in this expression represents as usual the average scattering by an isolated layer 
while the two following terms include the effects due to interferences. Sharp interference maxima 
may appear when a proper value Q“ approaches unity provided that the effect is not canceled by 
the vanishing of the R“* value in the numerator. 


(30) 


SECTION 5 


The results obtained in the previous section will be now applied to calculate interferences for 
mica type crystals with an arbitrary amount of randomness in the succession of the layers. It will be 
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assumed that each layer can occupy one of three positions with equal probability. The positions 
differ from each other in translations parallel to the plane of the layer by one-third of an identity 
distance (primitive translation). The distances between neighboring layers are identical so that all 
phase shifts between neighboring layers may be set equal to a definite quantity yg. The phase shifts 
due to the translations parallel to the planes of the layers can be incorporated in the layer form 


factors. All f values are equal to 3. 


The layers shall be assumed to be similar on their two faces. Then the most general expression for 


A, consistent with the values of f is 


3 — (a+b) a b 
Ai= a 3 —(a+c) c 
b c 3 — (b+) 
The quantities a, b, and c are the probabilities for finding a pair of layers of the kind 1, 2, ---, 2, 3 


or 3, 1, respectively. The probabilities that any layer shall be followed by its own kind can be derived 


from these numbers. 
The matrix Q; is given by 


3 — (a+b) a b 
i= a 3 —(a+c) c 3 exp (—i¢), 
b c 3— (b+c) 


which symbolizes that each element of the matrix is multiplied by 3 exp (—i@). The elements of the 


transformation matrix O may be obtained from the system of linear equations 
3 exp (—ig)[O¢(4—(a+b)) +a0¢? +60 ]=QMOeD, 
3 exp (—i¢g)[aO¢” +06 (3 —(a+b))+cO ]=QO), 
3 exp (—i¢g)[b0¢” +cO’® +068 (3 — (b+c)) ]=QO, 


where Q“) are proper values of the secular equation 


3[$—(a+5) exp (—i¢) ]—Q”, 3a exp (—i¢), 3b exp (—i¢) 
3a exp (—i¢), 3[3 —(a+c) exp (—iv)]—Q', 3c exp (—i¢) 
3b exp (—i¢), 3c exp (—i¢), 3[4 —(b+c) exp (—ig) ]— 


One finds the three proper values 
Q =exp (—i¢), 
QO =[1—3(a+b+c)+3M'] exp (—i¢), 
Q® =[1—3(a+b+c) —3M'] exp (—i¢), 
where M =[a?+6?+c?—ab—bc—ca | while for the matrix O one finds 
1//3 1/3 1/3 
O=|1//6[1—-A] 1/V/6Le+Ac*] 1/r/6[e*+Ac]}. 
1//6[1+A] 1//6Le—Ac*] 1/r/6[e*—Ae] 
Here the following abbreviations have been used 
e=exp (2771/3), «*=exp (—277/3), 
A= (ae+be* +c) M-}, 
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The matrix O-' may be obtained from O by interchanging rows and columns and taking conjugate 
complex values. 
The matrix F is simply 





; O O 
F=|O 3 O|. 
OO 3 


Since the three kinds of layers differ only in translations the absolute layers of the layer form 


factors are all equal and can be denoted by V. Corresponding to the translations by 3 primitive 


distance the layer form factors will differ by complex cube roots of unity and the V matrix can be 
written as 


! V2 eV? eV? 
| V=|eV? Ve eV?). 
eV? eV? Vp? 
Taking the diagonal elements of the matrix OV FO™ one obtains 
RY=0, R= RO =1y?, 


Thus the expression for the average intensity is 


QI [1—3(a+b+c)+3M'] exp (—i¢) [1—3(a+b+c)—3M!] exp (—i¢) ; 
—=2+ + +conjugate 
V2 1—[1-—3(a+b+c)+3M']exp(—ig) 1—[1—3(a+b+c)—3M!]exp (—ig) complex, 


(1-—D)(1+D-—F cos ¢) 
1+ F?+ D?—2F(1+D) cos 9+2D cos 2¢] 











Iu = r| (31) 


where the quantities F and D are defined by 


F=2-—6(a+b+0), 
0. D=1-—6(a+6+c)+27(ab+ac+bdc). 


For the case of complete order a=b=c=0, which leads to F=2, D=1. Then the numerator of 
Eq. (31) vanishes and there is no scattering except where g=n X2z in which case the denominator 
of Eq. (31) is also zero. By assuming small but not vanishing values for a, b, and c it can be shown 
that intensity maxima occur for g="X2z and these maxima go toward infinity when a, b, and c 
go toward zero. 

The case of complete disorder is obtained for a=b=c=}3. For this case the matrix A; shows that 
any sequence of two layers is equally probable. For this case F=D=0 and the intensity becomes 
equal to V? independently of ¢ as should be indeed expected from the discussion given in the beginning 
of Section 3. 





A third extreme case is obtained by setting average intensity becomes 
a=b=c=}. This is the maximum value that the 15+12 cos ¢ 
Ty —_ | | 





three quantities can attain if they are equal to 
each other. It physically corresponds to the situ- 33+40 cos +8 cos 29 

ation where a layer can be followed with equal This function has relatively broad maxima when 
probability by either of the other two layers but ¢g is an odd multiple of + and minima when 


not by itself. Then we have F= —1, D=}. The g=n X2z. The intensity in the minima is § of 
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the maxima. It is none too surprising that we 
now obtain minima where for an ordered lattice 
sharp maxima would occur since the scattering 
phases of neighboring layers are shifted 120° 
compared with the scattering of layers in an 
ordered structure. 

A simple example in which a, 0, and ¢ are not 
equal is furnished by a=%, b=c=0 which gives 
F=1, D=0. 


In = V?(1—cos ¢)/(2—2 cos ¢). 


This gives V?/2 whenever g¥#n X27. If p=n X24 
the expression becomes indeterminate. Assuming 
small but not vanishing values for } and c it can 
be shown that sharp intensity maxima occur for 
such values of ¢. The case discussed at present 
corresponds to a crystal broken up into two kinds 
of regions, first into a random mixture of the 
ist and 2nd kind of layer and second into an 
ordered structure of the 3rd kind of layer. The 
intensity function gives a structureless con- 
tinuum with sharp interference maxima super- 
imposed upon it. It must be noticed, however, 
that the sharp maxima are only partly due to the 
ordered portions of the crystals while they partly 
are caused by the fact that in the random parts 
of the crystal the average value of the scattering 
power is different from zero. 

It is of interest to notice that while the corre- 
lation probabilities are described by three 
parameters a, 6, and c, only two combinations of 
these parameters F and D enter into the intensity 
distribution. 

From the cases discussed above the effect of 
the two parameters F and D becomes apparent. 
According to their different values interference 
maxima may become sharper or broader and also 
the background may have various intensities. 
But as shown by the last example there is no 
unique relationship between breadth of maxima 
and intensity of the background. A background 
may appear even if the maxima are quite sharp. 


SECTION 6 


In the last example that is to be treated cor- 
relations are assumed between layers that are 
second neighbors so that the treatment of the 
previous section is not directly applicable. How- 
ever, by the artifice of considering pairs of layer 
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as forming, one layer of this case can be reduced 
to the one for which the general treatment has 
been given. 

A physical realization of the arrangement of 
layers that will be considered is a close packed 
structure of spheres in which the layers in 
question are the most densely packed arrange- 
ment of spheres. One of these, namely the cubic 
packed arrangement, can be represented by suc- 
cessive layers which are displaced with regard to 
each other by a translation perpendicular to the 
plane and by 3 of a certain primitive translation 
parallel to the plane. The three different positions 
obtainable by the latter translation will be desig- 
nated by I, II, and III and the cubic arrange- 


ment is characterized by the sequence ---I, II, 
III, I, I, III, 1--- or equally by ---TI, IIT, II, I, 
III, Il, I---. 


The second well-known close packing of 
spheres is the hexagonal close packed structure 
which can be described by either one of the three 
sequences 

---T, II, 1, 11,1, Il--- 
-+ TL, ITT, U1, 11, 11, I11--- 
-+ + TIT, I, 111, 1, II, 1--- 


In these structures a layer is always followed 
by a layer of a different kind and as far as inter- 
action between neighbors is concerned there is 
equal likelihood of a layer being followed by 
either one of the two remaining kinds of layers. 
The cubic and hexagonal structures differ in so 
far that in the cubic structure second neighboring 
layers must also be of a different kind while in 
the hexagonal arrangement second neighbors 
are of the same kind. Arrangements will now be 
discussed having certain degrees of randomness 
in which first neighbors are always different but 
second neighbors have certain probabilities to be 
alike or unlike.’ 

A second physical realization of the sequence 
of layers to be discussed is afforded by graphite 
crystals. The layers of carbon atoms have an 
arrangement as shown in Fig. 2. Successive 
layers are shifted parallel to the planes of layers 
by one-third of the primitive translation. In the 
graphite structure as described one finds an 


70. S. Edwards, H. Lipson, and A. J. C. Wilson, Nature 
148, 165 (August 9, 1941) have observed structures involv- 
ing mixtures of hexagonal and cubic close packing in 
metallic cobalt. 
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Fic. 2. Atomic arrangement in a carbon layer of graphite. 
Atoms in adjacent layers can either be above A or B, but 
only above one of them in an ordered crystal. 


ordered sequence in which every second layer is 
alike. Since in other samples of graphite com- 
pletely random shifts along the planes of the 
layers have been observed it is conceivable that 
crystals can be found in which the normal } 
primitive translation between layers is strictly 
preserved while an arbitrary correlation exists 
between translations of second layers. 

In order to reduce formulas to one with cor- 
relations between neighboring layers pairs of 
layers are considered as units. Six pairs can exist 


(1, TI), MDC, Hd, DAM, I, 1D 


each of which occurs with the frequency }. First 
consider by what kinds of layers (I, II) can be 
followed. According to the above statements the 
sequences (I, II), (II, III) and (I, II)(II, I) 
cannot occur. It is assumed that the frequencies 
with which the combinations (I, II)(I, II), 
(1, ID(1, IT), and (1, ID)(1I1, 1) occur are in 





(II, Z)(11, IT) (1, 111) 


(II, J), M1), 11) 


(II, 7) (1H, 11) (1, 11) 


(II, Z)(I1I, Hd, 11) 
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the ratio 1:x: y. The combination (I, I1)(III, IT) 
has the same frequency as (I, II)(I, III). In fact 
the probabilities remain unchanged if the desig- 
nations I, II, and III are permuted in any way 
or if sequences are read in the opposite directions, 
since the original restrictions relate only to the 
questions of similarity or dissimilarity of first and 
second neighbors. Now the sequence (I, IT)(1, IIT) 
can be turned into the sequence (I, IT) (III, IT) 
by substituting I for III, II for I, III for II and 
reading in the reverse direction. By similar 
methods one can obtain all other sequences from 
those mentioned above. 

A connection between the relative proabilities 
x and y can be established if it is recalled that 
the grouping of the layers I, II, and III into 
pairs was arbitrary. Consider the question of the 
probability with which a layer I will be followed 
by the same kind of layer in the third place. 
There are two ways in which this can happen; 
namely, by the arrangement (I, II)(III, I) or by 
(I, I11)(II, 1). Assuming that the first layer is of 
the kind I the probability that the first pair of 
lavers is (I, II) will be 3. This pair can be fol- 
lowed by the combination (I, II), (1, III) (III, I1) 
X (III, 1) with the relative probabilities 1, x, x, 
and y. Therefore the probability for the sequence 
(I, I1)(1II, I) is 3Ly/(1+2x+y) ] supposing that 
the first laver was of the kind I. By a similar 
argument it can be shown that the sequence 
(1, I11)(11, 1) has the same probability so that 
the total probability of I being followed in the 
third place by I is y/(1+2x+y). 

The same probability may be calculated also 
in a somewhat more complicated manner by 
adding the probabilities of occurrences for the 
following eight combinations 








x x 
1H2x+y 142x453 
x Vy 
142xr+y 142x+y 
x x 
1+2x+y¥ 14+2x+y 
y y 
142r+y 142e+y 
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and four more combinations in which layers two and three are systematically replaced by each other. 
The probabilities of these combinations are indicated after the symbols assuming that the italicized 
layer is of the kind I. These four probabilities have to be added and multiplied by two to take into 
account combinations in which II and III are interchanged. One obtains [(x+y)/(1+2x+y) ? for 
the same probability for which y/(1+2x+~) was obtained above. These two quantities can be equal 
only if y=’. Thus all probabilities can be reduced to the relative probabilities 1 : x for the sequences 
(I, 11) (1, IT) and (I, IND (1, II). 

Taking into account the relative probability of sequences and the frequencies for individual pairs 
one obtains for the A matrix 


I,tl or $@lt U,t U1, I,1ll 








1,1} 1 0 x 0 x x 
11, 11] x2 1 0 x 0 Xx 
III, 1] 0 x 1 Xx Xx 0 1 
~ ILI 0 x Xx 1 X 0 lox+1)? 
Ill, | x 0 Xx 0 1 x? 
I,1t| xX x 0 x? 0 0 


where it is indicated that each element of the matrix is to be multiplied by 1/6(X +1)’. In the above 
form it has been indicated to which pair of layers each row and column corresponds. It will be ob- 
served that A is not symmetrical. 

The phase shifts between consecutive pairs of layers may be written as the same number ¢ for any 
sequence. The effect of shifts parallel to the layers may be included in the complex layer form factors 
V. The matrix Q; then becomes 








1 oo x2 0 xX Xx 

xX? 1 0X 0 X 

exp(—ig)|} 0 X? 1 X xX O 
atx |0 xX xX 1 xX? 0 
Xx 0 X 01 X 

x xX 0 X: 0 1 


where exp (—tv)/(1+X)? written before the matrix indicates that each element of the matrix contains 
this as a factor. 
This matrix can be diagonalized by means of the matrix 


; £ 8 1 1 1 
11 1 —1 —1 —1 
1ji ¢ & [at] [at]e [at let 
J/6l1 & ¢€ [at]}* [at]*e* [at]*e 
le & [a] [o]Je  [al]e* 
1 ¢€ [a }* [an ]*e* [am ]*e 








and by the reciprocal matrix O-, 












































ains 





PARTIALLY ORDERED LATTICES 165 


The symbols ¢ and e* have the same meaning as in the first part of the section and the abbreviations 
x+(4—3x?)} 
2 
x— (4—3x?)! 


eS + ex. 
2 


at +ex, 


Their conjugate complex values are indicated by asterisks. With the transformation OQ,0~' the 
proper values of Q; are 
Q =exp (—7¢), 


x—1\? 
g0=(5=2)' exp io 
x 


1 — [Lat |x + ex? 














0% =Q® = ecm exp (—1¢), 
: x a 
1—[am Jv+ex*) 

Q® =Q® = ; - oe ~ exp (—i¢). 
a x}° a 


It will be noticed that the matrix Q; could be diagonalized in spite of the fact that A is not sym- 
metrical and all the proper values of Q are not different. Actually the fulfillment of one of these con- 
ditions as given in Section 4 is a sufficient but not necessary condition for transforming Q; to a 
diagonal matrix. 

The matrix F is a diagonal matrix with all diagonal elements equal to §. The matrix V may be 
written after taking into account the translation by 3} and 3? the primitive translations in layers of 
type II and III 

| 7+ eT+ eF*+ T eT eT | 
e*T+ Tt ¢eT+ &T T eT 
T+ T+ T+ GT @T T | 

T eT eT T- eT lk 

e*T T a £a°”lUC<“Cda;llCO|SS 


eI TT T . @ FF" 


V=2|V/2 





where the abbreviations 
T= —}+ cos 9/2, 


T+=1+cos (¢/2+27/3), 
T-=1+cos (¢/2—27/3) 
have been used. The quantity | V|* is the square of the absolute value of the layer form factor for a 


layer of the type I, II, or III. 
Calculating the diagonal elements of the matrix 


R=OVFO" 
one obtains 
ROD = RO = RO = RO =, 


[a~]T*++[at]7-+27 
[at]—La-] 

[at ]T+—[a-]T--—2T 

[at]—-[o-] 


Re = — 





Ro = 
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Fic. 3. Scattering from mixed cubic (x= ©) and hexagonal (x=0) close packed structures. The solid curve gives the 


value of expression (34) for x =3.0 and the dotted curve for x=0.3. 


Substituting these quantities into Eq. (30) the average intensity per laver is found to be 
In (—Ca-JT++[ot]T-+27)(1—Lat Jeter’) exp (—ig) 
Vi? (Cat ]—[La-)[(1+x)?— (1 — [at ]x-+ex’) exp (—ig)] 

(Cat ]T+—[a-]T-—2T)(1—[a- ]x+ex*) exp (—i¢) 

* (Cat ]—[a-])C[(1+x)?—(1—[Lam Jx+ex?) exp (—7¢) ] 








+conjugate complex +(7++77-). (32) 
The last term in Eq. (32) gives the average scattering of an isolated layer. Equation (32) can be 
expressed in the somewhat simpler form 
Ty _ (1+)? [ —[a-~]T++[at]T7-+2T 
[V2 Cat ]—Lo-J(1+2)*—[1 [ar }rt er?) exp (—ig) 
+[at]T+-—[La-]T--—2T 
(1+x)?—[1—[aq ]jx+ex?] exp (—7¢) 


Substituting into Eq. (33) the values of [at], [a~], [7+], and T~ one finally obtains 








|-+eonjusate complex —(2—cos ¢/2). (33) 


In (1+)? p (4 —3x*)*(2 —cos ¢/2)+3ix sin ¢/2—2+4 cos ¢/2 
Vi? (4—3x2) — 2(1-+x)?—(2—x* —x(4—3x%)!) exp (—iy) 





(4—3x?)!(2—cos ¢/2) —3ix sin ¢/2+2—4 cos ¢/2 _ 
ns = [+eonjugate complex —(2—cos g/2). (34) 





2(1+x)?—(2—x?+x(4—3x?)!) exp (—7¢) 


The limiting case of hexagonal close packing is obtained for x =0. Substitution in (34) shows that 
T,, indeed vanishes for x=0 and exp (ig) #1. Now exp (ig) =1 is just the interference condition for 
a hexagonal lattice since g was the phase shift between second neighbors and this is the identity 
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distance in the hexagonal packing. For exp (tig) =1 and x=0 Eq. (34) becomes indeterminate. By 
carrying out the process x—0 one obtains intensity maxima in the correct position that approach 
infinity as 1/x. 

The other limiting condition x= ~ corresponds to a cubic close packed arrangement. In this case 
Kq. (34) vanishes except when ¢ is an integral multiple of 47 3 but not of 27. Indeed if second 
neighbors have such phase shifts the phase shift of third neighbors becomes a multiple of 27. For 
cubic close packing the distance of third neighbors is the primitive distance. If y is a multiple of 
47/3 but not of 27 and if x <« Eq. (34) approaches ~ proportional to x and we again obtain sharp 
interference maxima. These maxima may be interpreted as belonging to a mixture of two cubic close 
packed crystals; one with the sequence of layers ---I IL HI 111 III--- gives interference maxima 
at g= —4nr/3+nX2r (n is + — or 0); the sequence of layers ---I II] IT I III II--- gives maxima 
at g=4n/3+nX2r. 

It can be seen from the form of Eq. (34) that for values of x between zero and infinity no sharp 
maxima occur and indeed these intermediate values correspond to imperfectly ordered arrangements. 
Values of Eq. (34) calculated for x=0.3 and x=3.0 are shown in Fig. 3. Positions of interference 
maxima for the ordered systems having x=0 and x= ~ are also indicated. 
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Supersonic Measurement of the Heat Capacity of Propylene 


D. TELFAIR* 
Department of Physics, The Pennsylvania State College, State College, Pennsylvania 
(Received November 22, 1941) 


An apparatus designed for (1) a general study of acoustical dispersion in gases, and (2) meas- 
urement of gaseous heat capacities has been used to measure the heat capacity of propylene 
between 273 and 490°K. The general methods used are outlined and discussed. The measured 
heat capacity values agree weil with those calculated from the frequency assignment as given 
by Wilson and Wells, and an assumed potential restricting rotation of about 2000 cal./mole, 
but deviate from the values calculated from Pitzer’s frequency assignment and an 800-cal./mole 
potential. The measurements on propylene also have the value of serving as a test of the 
usefulness and practicality of the supersonic method of heat capacity determinations. 


INTRODUCTION ever, made it apparent that sonic methods might 
yield misleading results, since the presence of 
lagging vibrational states is accompanied by an 
abnormally high velocity unless the frequency 
used is low enough to allow all modes of storing 
energy to stay in equilibrium with the fluctuating 
pressure and temperature which accompany the 
passage of a sound wave through a gas. 


ARIOUS workers have used sound velocity 
measurements as a means of determining 
gaseous heat capacities.'~7 For some time this 
was considered as one of the most precise and 
reliable methods. The discovery and subsequent 
investigation of intramolecular dispersion, how- 





_*Now with the Plastics Division of the Monsanto The development of the Pierce acoustic 
Chemical Company at Springfield, Massachusetts. . 8 : 
' Kundt, Pogg. Ann. 135, 337 (1868). interferometer’ put into the hands of the 
*Wullner, Ann. d. Physik 4, 321 (1878). experimenter an instrument capable of very 
* Capstick, Phil. Trans. Roy. Soc., 185, 1 (1895). : “ae : = : ’ 
‘E. H. Stevens, Ann. d. Physik 7, 285 (1902). precise velocity determinations, but only a few 
* Kalahne, Ann. d. Physik 11, 225 (1903). attempts have been made to use it for heat 
12) Campbell, and Parker, Proc. Roy. Soc. 100, 1 


7 J. R. Partington and W. C. Shilling, Phil. Mag. 3, 273 8G. W. Pierce, Proc. Am. Acad. Arts and Sci. 60, 271 
(1927); Phil. Mag. 6, 920 (1928). (1925). 
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capacity measurements*" principally for the 
reason outlined above. 

However, it seems probable that the sonic 
method of heat capacity determinations can 
still be used in many cases, especially when 
viewed in the light of recent work by Eucken and 
his co-workers.’ '!5 The latter have pointed out 
that for all polyatomic gases examined so far, 
one obtains a single smooth dispersion curve 
similar to that for the simple case of an idealized 
gas with only one vibrational quantum state 
above the ground state. Thus, gases behave as 
if all of the vibrational quantum transitions have 
the same (or nearly the same) relaxation time. 
By locating the corresponding single dispersive 
region with preliminary velocity and absorption 
measurements, one may frequently choose the 
conditions of measurement (temperature, pres- 
sure, and frequency), in such a way as to be on 
the normal, low frequency, low velocity side of 
the dispersive region, and consequently measure 
a velocity corresponding to the true, or total 
heat capacity. The application of this method 
to gaseous propylene is described below. 


METHODS USED 


The apparatus used will be described in detail 
in another journal.} Briefly, it consists of a 
Pierce acoustic interferometer which, in general 
design and purpose is similar to the one described 
by Eucken and Niimann.'* The interferometer 
chamber, which may be evacuated, is operated 
in conjunction with a temperature control 
system, a vacuum system, and a potentiometer 
for temperature measurement. The half-wave- 
length spacings are accurately determined with 
a calibrated micrometer screw, and the frequency 
of the quartz crystal source is determined with 
high precision by measuring the audible beat 


° 0) C. Shilling and J. R. Partington, Phil. Mag. 9, 1020 
(1930). 

10 A. Eucken and R. Becker, Zeits. f. physik. Chemie 
B27, 235 (1934). 

1 R. Schultze, Ann. d. Physik 34, 41 (1939). 

2S. K. K. Jatkar, Ind. J. Phys. 13, 445 (1939). 

18 Van Lammeran, Physica 2, 833 (1935). 

14 A. Eucken and S. Aybar, Zeits. f. physik. Chemie B46, 
195 (1940). 

1 K. Schafer, Zeits. f. physik. Chemie B46, 212 (1940). 

+ D. Telfair and W. H. Pielemeier, Rev. Sci. Inst., in 
press. 

16 A. Eucken and E. Niimann, Zeits. f. physik. Chemie 
B36, 163 (1937). 
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note obtained with an ordinary broadcast 
receiver when the radiofrequency signal from 
the crystal is picked up together with a broadcast 
carrier signal of frequency close to that of the 
crystal. The velocity is then obtained from the 
product of the wave-length and the frequency. 

The absorption is calculated from plate 
current readings taken simultaneously with the 
micrometer screw readings which are used in 
calculating the velocity. The peak and trough 
values of the plate current are recorded for each 
setting of the micrometer screw, and the ab- 
sorption is calculated from the logarithmic 
decrement of these peak-to-trough differences of 
the plate current, as the reflector is moved away 
from the crystal. This method has been justified 
experimentally by Pielemeier'? and van Itter- 
beek'® and has been given a much more satis- 
factory theoretical basis by Hardy,'® who has 
demonstrated that it gives a good approximation 
to the true absorption, provided a small ampli- 
tude source and path differences far from the 
source are used. These conditions were satisfied 
for the measurements reported below. 

The velocity of propagation of plane sound 
waves in any fluid is given with sufficient 
exactness for all practical cases by the equation 


E; C. v Op 
V2=y—= -=-(—) : (1) 
p - Cy, p\ dv/ 7 


where 7 =the heat capacity ratio, C,/C., 
E;=the isothermal elasticity, 
p=the density, 
p=the pressure, 
v=the volume. 


The equation of state may be used (a) to 
evaluate the isothermal elasticity, —v(dp/dv)r, 
and (b) to express C, in terms of C,. Thus, for 
a perfect gas, a useful form of Eq. (1) becomes 


V2?=(RT/M)(1+R/C,), (2) 


where R=the gas constant per mole, 
M=molecular weight, 
T =the absolute temperature. 


17 W. H. Pielemeier, Phys. Rev. 34, 1184 (1929). I 

18 A. van Itterbeek and P. Mariéns, Physica 4, 207 
(1937). 

19H. C. Hardy, Doctor’s Thesis, Pennsylvania State 
College (1941) (to be published soon). 
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Fic. 1. Measured and calculated values of the constant 
volume molal heat capacity of propylene. 


This same useful form may be retained for real 
gases. By assuming an equation of state of the 
form 


pu= RT +8/0+6/v?+e€/v+---, (3) 


Hardy” has carried through a rigorous derivation 
of the equation 


V?=(RT/M)(P+QR/C,). (4) 
The quantities 


P=1+28/RTv+36/(RTv*)+:-- (5) 


and 
1 0B 1 06 2 
o-[14+— (—) i —) + | (6) 
Rv\0T/, Rv? \dT/, 
are dimensionless and approach unity as the gas 
approaches perfection. Thus, if the velocity has 
been measured, the molal heat capacity at 


constant volume may be obtained from Eq. (4). 
By solving for C, we obtain 


OR 
" (V2M/RT—P) 


- 





(7) 


Finally, using the value of C, calculated by 
means of Eq. (7), we obtain C,°, the constant 
volume molal heat capacity for an infinitely 
dilute gas, from the familiar thermodynamic 


expression 
v 0*p 
C0=C.— f r(—) dy, (8) 
Ps éT*/, 


*°H. C. Hardy, Paper No. 1, Rochester Meeting of the 
Acoustical Society of America, May 5-7, 1941. 
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which, using (3) as the equation of state, becomes 


TT s0°8 1 £076 
ca=c.+—|(—) bine (—) +] (9) 
vL\oT?/, 2 N\oaT?/, 


If absorption measurements are made simul- 
taneously with velocity measurements, a con- 
sideration of the way in which the absorption 
changes with frequency and/or gas density will 
show on which side of the dispersive region these 
measurements lie. Also, the numerical value of 
the absorption allows one to estimate quantita- 
tively how far the measurements are removed 
from the dispersive region. 


THE HEAT CAPACITY OF PROPYLENE 


The propylene used in these measurements 
was furnished by the Pennsylvania State College 
Cryogenic Laboratory. Vapor pressure readings 
taken during its purification indicated that the 
impurity was not greater than 0.05 percent. 

Preliminary measurements of velocity and 
absorption taken at various frequencies, pres- 
sures, and temperatures show that even for 0°C 
and 1.2 megacycles, propylene is well on the 
normal low frequency, low velocity side of the 
dispersive region. The absorption, measured as 
described above, is very nearly as low for 
propylene at all the temperatures at which 
measurements were taken, as it is at room 
temperature for air, or nitrogen, which have 
very little intramolecular absorption. Conse- 
quently, propylene, under these conditions, must 
be far removed from any dispersive regions. 
By making use of (1) the measured absorption, 
(2) the fact that the absorption increases with 
frequency, and (3) the expression which gives y, 
the absorption per wave-length, as a function of 
the frequency,”' it is possible to calculate an 
approximate value for the frequency at which 
maximum absorption would occur in propylene. 
If based on the reasonable assumption that all 
modes of vibration in propylene have the same 
relaxation time, such a calculation indicates that 
the frequency for maximum absorption would 
occur between 80 and 100 megacycles. This 
corresponds to a relaxation time of approxi- 
mately 10-8 second. 


*1W. H. Pielemeier, H. L. Saxton and D. Telfair, J. 
Chem. Phys. 8, 106 (1940). 











170 


From these considerations we may conclude 
that either (a) all of the quantum vibrational 
states of propylene, under the conditions of our 
measurements, have such a short mean quantum 
lifetime that they are able to stay in equilibrium 
with the rapid pressure and temperature fluctua- 
tions of the sound wave, or (b) only states 
corresponding to high wave numbers, and 
therefore contributing only slightly to the total 
heat capacity, have “dropped out’ of the 
acoustical cycle at much lower frequencies. (If 
states of low wave numbers had ‘dropped out”’ 
at lower frequencies, the absorption would, of 
necessity, decrease with an increase in log (f/p), 
whereas it is observed to increase.) 

As evidence that (a) actually obtains rather 
than (b), we present the close agreement between 
the values of C,° (the constant volume molal 
heat capacity for infinite dilution) as calculated 
from velocity measurements with those deter- 
mined by Kistiakowsky and Rice” by an 
independent method (see Fig. 1). 

Berthelot’s equation of state was first used for 
the C,° calculations but was found to be entirely 
inadequate. Fortunately, the second virial coeff- 
cient 6 has been determined for propylene 
experimentally. We have used the value given 
by Roper* in order to calculate C,° for each of 
the seven temperatures indicated in Fig. 1. Since 
small ranges of pressure were used (from } to 2 
atmos.), the second virial coefficient only was 
considered as sufficient for our calculations. 

In order to present the results in a more 
convenient form, the C,° values as read from a 
larger scale plot of Fig. 1° for temperature 
intervals of ten degrees lying between 270 and 
510°C have been listed in Table I. uf 


THE RESTRICTED ROTATION POTENTIAL AND 
FREQUENCY ASSIGNMENT IN PROPYLENE. . 


As recently pointed out,™:a_ potential’ re- 
stricting rotation of the methyl group in pro- 
pvlene of about 800 cal./mole, and the funda- 
mental vibrational frequencies as assigned by 
Pitzer®® do not disagree with the measured heat 

22 G. B. Kistiakowsky and W. W. Rice, J. Chem. Phys. 
8, 610 (1940). 

23 E. E. Roper, J. Phys. Chem. 44, 835 (1940). 

* TD. Telfair and W. H. Pielemeier, J. Chem. Phys. 9, 571 


(1941). 
% KX. S. Pitzer, J. Chem. Phys. 5, 473 (1937). 
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capacity of propylene between 273 and 367°I< 
as reported by Kistiakowsky and Rice.” How- 
ever, the set of frequencies recently assigned by 
Wilson and Wells,?® and a potential restricting 
rotation of about 2000 cal./mole are also in 
agreement with the measured values of the heat 
capacity. By extending the heat capacity 
measurements to higher temperatures we have 
made possible a choice between the two assign- 
ments. As shown in Fig. 1 the measured heat 
capacities at the higher temperatures (between 
360 and 500°K) are still in excellent agreement 
with values calculated from the Wilson-Wells 
assignment and a 2000-cal./mole potential, but 
values calculated from Pitzer’s assignment and 
an 800-cal./mole potential deviate by consider- 
ably more than the limits of accuracy of the 
measurements. In these calculations we have 
used Pitzer’s tables’? for the heat capacity of a 
restricted rotator. Pitzer and Gwinn** have 
recently shown that these tables are applicable 
to a molecule like propylene, which consists 
essentially of a rigid frame with an attached 
symmetrical top whose rotation is hindered. 


TABLE I. Constant volume molal heat capacities of pro- 
pvlene for an infinitely dilute gas. 











C, | CY 
Temperature cal/mole Temperature cal/mole 
270°K 12.30 400°K 17.17 
280 12.70 410 17.53 
290 13.09 420 17.89 
300 13.48 430 18.24 
310 13.86 440 18.59 
320 14.24 450 18.94 
330 14.62 | 460 19.29 
340 14.99 | 470 19.63 © 
350 15.36 480 19.98 
360 15.73 | 490 20.32 
370 16.09 500 20.66 
380 16.45 510 21.00 
390 16.81 








. Attention is called to the fact that the slope 
of the Wilson-Wells curve is everywhere in 
agreement with the slope for the measured values 
of the C,° curve, whereas the Pitzer curve is not. 
Thus, Pitzer’s frequency assignment with a 
potential barrier lower than 800 cal./mole might 
26 E. B. Wilson and A. J. Wells, J. Chem. Phys. 9, 319 
(1941). 

27K. S. Pitzer, J. Chem. Phys. 5, 469 (1937). . 

28K. S. Pitzer and W. D. Gwinn, J. Chem. Phys. 9 455 
(1941). 
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give better agreement in the higher temperature 
region, but only at the expense of the present 
good agreement in the lower temperature region 
of the curve for measured C,° values. 

It should also be mentioned that the entropy 
of propylene calculated on the basis of the 
Wilson-Wells frequency assignment and a 2000- 
cal./mole potential with the 
measured entropy provided one assumes no 
entropy of mixing in the solid. On the other 
hand, if an entropy of mixing is assumed, the 
entropy calculated on Pitzer’s basis agrees with 
the measured value. However, judging from the 
results of a study of cyclopentane by Aston and 
his co-workers,”® an entropy of mixing cannot be 
expected in solid propylene. 


barrier agrees 


ESTIMATED PRECISION OF MEASUREMENTS 


If the wave-length frequency product Af is 
substituted for the velocity V in Eq. (7) one 
obtains 


OR 
‘(Af)2M/RT—P}_ 


C,= (10) 


The error of measurement present in C, depends 
on the error of measurement in the three meas- 
ured quantities A, f, and 7. One disadvantage 
of the sonic method of heat capacity measure- 
ments lies in the fact that a rather large propa- 
gation of errors is involved in the calculation of 
C,. The quantity (Af)?//RT is approximately 


* J. G. Aston, J. Am. Chem. Soc. 63, 2029 (1941). 
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equal to the heat capacity ratio, C,,C,. Thus 
(since P is close to unity) the denominator of 
Eq. (10) will be from one-fourth to one-twelfth 
as large as (Af)?//RT itself, and consequently 
the relative error entering into Eq. (10) will be 
from four to twelve times as large as the relative 
error in (Af)?M/RT. For the propylene measure- 
ments, each of the quantities f?, \*, and 7 were 
measured with a precision of 1 part in 2000, or 
better, and consequently the resulting C, values, 
on the average, would be in error by about 
0.6 percent because of errors in A, f, and 7. If 
additional precautions are taken in general 
technique, greater precision than this 
should be possible. 

Another quantity which may introduce error 
into the C, determinations is uncertainty in the 
molecular weight M. As already stated, a 
specially purified sample of propylene was used 
for these measurements, and care was taken to 
retain its purity while measurements were being 
taken. 

The 1937 international atomic weights (C 
=12.010 and H=1.0078) have been used. The 
ice-point has been taken as 273.16°K, and for 
the gas constant 8.3415 X10" ergs/deg. mole, or 
1.9869 cal./deg./mole has been used. 

The author wishes to express his appreciation 
of the helpful criticism and guidance in this work 
given by Professor W. H. Pielemeier and 
Professor J. G. Aston, and of the stimulation 
and encouragement of his colleague, Dr. H. C. 


Hardy. 
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Raman displacements, estimated intensities and depolarization factors have been obtained for 
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CH;CH2C =CH, CH;(CH2)2C =CH, CH;3(CH2)s3C =CH, CH;CH2C =CCH2CH:, CH3(CH2)2C 
=C(CH2)2CH;, and CH;(CH2)4C=CClI. The data obtained in the 2200 cm™ region were 
[2068(1)P], 2118(10)0.5, 2158(3)0.5; [2064(0)], 2118(10)0.5; 2071(0), 2094(0), 2118(10)0.4; 
2200(1)P, 2231(9)0.4, 2247(8)0.4, 2301(6)0.4; 2234(10)0.5, 2292(6)0.5; 2148(1)P, 2194(1)P, 
and 2243(10)0.4, respectively. Unusual features are the appearance of the second strong line 
at 2158 cm™ in the spectrum of 1-butyne and the absence of the line near 2245 cm in the 
spectrum of 4-octyne. Introduction of the Cl atom in the last compound radically modifies the 


spectrum in this region. 








INTRODUCTION 


ISUBSTITUTED acetylenes have more 
lines in the 2200-cm— region of the Raman 
spectrum than monosubstituted acetylenes. The 
present series of investigations! has been directed 
toward a better understanding of the reasons for 
this phenomenon. In particular, it was desired to 
observe the changes produced in this region when 
various substituents replaced one or both of the 
hydrogen atoms in the acetylene molecule. 

In I and II, studies were made of a series of 
disubstituted acetylenes of the type CsH;C =CR. 
It was found that the spectrum in the 2200 cm 
region is quite sensitive to changes in the sub- 
stituent R, the spectrum being dependent not 
only upon the atom or group immediately at- 
tached to the acetylenic carbon, but also upon 
the atom or group once or twice removed from it. 
In general, however, at least two strong lines 
with approximately equal depolarization factors 
were observed near 2215 and 2240 cm—. 

Later, an infra-red investigation? of CsHsC 
=CCH; was made in an attempt to test the 
Badger* hypothesis that Fermi resonance inter- 
action was responsible for some of the additional 
lines in disubstituted acetylenes. While the fre- 
quency conditions for resonance did not seem to 


‘I, J. Am. Chem. Soc. 60, 2664 (1938); II, ibid. 61, 3546 
(1939); III, bid. 62, 3185 (1940); IV, J. Chem. Phys. 9, 


Chem. 


390 (1941); V, J. Am. Chem. Soc. 63, 1 18 (1941). 
2M. Murray and Forrest F. Cleveland, J. 
Phys. 8, 133 (1940). 
7R. M. Badger, J. Chem, Phys. 5, 178 (1937). 
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be very satisfactorily fulfilled, it was difficult to 
draw definite conclusions because of the uncer- 
tainty in regard to the effect of the anharmonicity 
upon the overtone frequencies. 

In III, CsH;C =CI and eight substituted acety- 
lenes not containing a phenyl group were studied. 
For the monosubstituted acetylenes, it was found 
that the substitution of a tertiary carbon atom 
caused a slight decrease in the 2118 cm~ fre- 
quency and that in the case of CH;C(CH;)(OH)¢ 
=CH two strong lines, instead of the usual one 
for monosubstituted acetylenes, were observed 
at 2089(3) and 2112(10). Single lines were ob- 
served for two of the disubstituted acetylenes, 
but this is probably not significant since the line 
in one case was so broad as to suggest an un- 
resolved doublet and in both cases the spectra 
were incomplete and a weak line could have gone 
unobserved. 

In IV, it was shown that some of the weaker 
lines in the spectra of both monosubstituted and 
disubstituted acetylenes, having displacements 
smaller than that of the most intense line, may 
be due to isotopic molecules having C® at one of 
the triple-bond positions, as originally suggested 
by Glockler and Renfrew for 2-butyne.* 

An interesting case® is that of CH 3(CH2)sC 
=C(CH.)4,CH; in which the substitution of an 
oxygen atom for two hydrogen atoms to give 


958) Glockler and M. M. Renfrew, J. Chem. Phys. 6, 408 
(1 3 

5M. J. Murray and Forrest F. Cleveland, J. Am. Chem. 
Soc. 63, 1363 (1941). 
















RAMAN SPECTRA OF ACETYLENES 


CH3(CH2)4C =CCO(CH2)3CH; changed the spec- 
trum from 2231(7), 2248(2), 2294(4) to a single 
line at 2212(7). 

In V, several alkyl acetylenes were investi- 
gated. The two monosubstituted ones gave the 
strong line at 2118 cm~ and also weak lines near 
2064 and 2095 cm~', which may be due to the 
C® isotopes. All of the disubstituted acetylenes 
gave a strong line near 2235, a slightly weaker 
one near 2295 and all except 9-octadecyne, for 
which the spectrum was incomplete, gave a weak 
line near 2250 cm-!. 

The present paper gives the results of an 
investigation of several of the lower members of 
the series of alkyl acetylenes. It was hoped that 
this investigation might provide some clue as to 
the origin of the line at 2250 cm~. 


EXPERIMENTAL 


The apparatus and experimental technique 
have been described® previously. The spectra 
were produced by use of Hg 4358A. Details in 
regard to the preparation and physical properties 
of the compounds are as follows: 


1-Butyne, CH;CH.C =CH, prepared from so- 


dium acetylide and ethyl bromide according to 
the method of Vaughn, Hennion, Vogt, and 
Nieuwland,’ b.p. 18—21°C. 

1-Pentyne, CH;(CH2)2C =CH, prepared analo- 
gously, b.p. 41—43°. 

1-Hexyne,*CH;(CH2)3C =CH, prepared analo- 
gously, b.p. 72°. 

3-Hexyne, CH;CH:C=CCH.CHs, prepared 
from sodium acetylide, sodamide and ethyl 
bromide, according to the method of Bried and 
Hennion,? b.p. 80.5-81.0°. . 

4-Octyne, CH3(CHe)eC =C(CH2)sCHs, pre- 
pared analogously, b.p. 130.0-130.5°. 

1-Chloro-1-heptyne, CH;(CH:2),C=CCl, from 
1-heptyne to CsHi.C =CNa by sodamide in 
liquid ammonia, thence, after removal of am- 
monia, to C;H;,;C=CCI by action of CH3C.H, 


° J. Chem. Phys. 7, 396 (1939); ibid. 8, 153 (1940); ibid. 
10, 18 (1942). 

7*T. H. Vaughn, G. F. Hennion, R. R. Vogt, and J. A. 
Nieuwland, J. Org. Chem. 2, 1 (1937). 

* Synthesized by Mr. R. E. Dineen and spectrograms 
measured by Mr. Robert D. Sieg. 

°E. A. Bried and G. F. Hennion, J. Am. Chem. Soc. 59, 
1310 (1937). 
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SO2CI in ether, b.p. 58.0-58.4° at 35 mm pres- 


sure. 
RESULTS 


The results are given in Table I. The only 
compounds for which data seem to have been 
reported previously are 1-pentyne,'® 1-hexyne!® ! 
and 3-hexyne." Seven new lines were observed 
for 1-pentyne, and the 1446 line was resolved 
into the doublet 1432, 1454. The mean deviation 
of the present and former frequencies was 4 cm~', 
the maximum deviation being 12 cm~' (for the 
758 and 1044 lines). For 1-hexyne, eight new 
lines were found, the doublet 801, 814 appeared 
unresolved at 810 and the lines at 431 and 1024 
were not found. The mean and maximum devia- 
tions were 2 cm~'! and 9 cm~! (for the 2908 line), 
respectively. And finally, eleven new lines were 
observed for 3-hexyne, the 492 and 1447 lines 
were resolved into the doublets 484, 499 and 
1439, 1460, the doublet 196, 211 appeared unre- 
solved at 206 cm and the 321 and 1144 lines 
were not found. The mean and maximum devia- 
tions were 1 cm~! and 3 cm~', respectively. 
Depolarization factors do not seem to have been 
reported heretofore for any of the present com- 
pounds. 


DISCUSSION OF RESULTS 
Frequencies Near 2200 cm™ 


Monosubstituted Acetylenes 


The Raman data obtained in this laboratory 
for a number of monosubstituted acetylenes are 
collected in Table II; 1-propyne is also included 
in the table for completeness. The weak lines 
near 2066 and 2095 may be ascribed to the 
triple-bond frequencies of the two possible C™ 
isotopes (see IV). The strong line at 2124 cm™! 
for 1-propyne drops to 2118 for 1-butyne and 
retains this value for higher members of the 
series, as well as for many other monoalkyl- 
acetylenes. For 1-butyne, an additional line of 
medium intensity was observed at 2158 cm™. 
This line may be due to Fermi resonance inter- 
action between the triple bond fundamental and 


10 R. Courtel, Dipléme d’études supérieures, Paris, 1932. 

u . Bourguel and P. Daure, Bull. Soc. Chim. 47, 1349 
(1930). 

2H. Van Risseghem, Bull. Soc. Chim. Belg. 47, 271 
(1938). 
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TABLE I. Raman spectra of some acetylenes. 


1-Butyne ntyne 1-Ilexyne 


I p Av I 


5b 


1094 


1228 
1261 


1300 
1327 


1432 
1454 
[2064 


2118 


P 274 : 2737 
0.4 ¥s 2839 
0.4 , 2869 
0.2 ‘ od 2908 
0.2 : 2933 
0.8 s , 2965 
0.6 3304 3303 


2851 
2885 
2919 
2941 
2986 
3305 


NUON Wwe 
co 


1-Chloro-1- 


3-Hexyne 4-Octyne Heptyne 


Av I 


0.6 


0.8 
0.9 
0.8 


[0.9] 


1020 
1037 
1069 0.7 
1095 4 0.9 

0.6 


{0.9} 


1227 
1260 
1274 
1291 


1327 
1350 


1434 


P 
0.5 
0.4 
0.5 


P 27 P 


03 j 03 
03 ‘ 0.5 0.3 
0.1 K 0.4 K 0.2 
0.7 297 0.9 2965 3 0.9 








the first overtone of the strong, polarized Raman 
line observed at 1068 cm~'!. However, several 
combinations of the observed Raman frequencies 
fall in this region and it may be that one of these 
combination frequencies is involved in the reso- 
nance interaction that accounts for the additional 
line. Another possibility is that an overtone of a 
vibration (or a combination of two vibrations) 





Av=Raman displacement in cm~!, J =estimated intensity, p =depolarization factor, tindicates observation both as Stokes and anti-Stokes 
lines, brackets enclose data in regard to which there is some question, b =broad, vb =very broad, P means that p was less than 6/7 but that it 
was not possible to make a quantitative measurement, the brace joins lines that were unresolved on the polarization spectrogram, and bold type 
indicates previously observed frequencies. 


that is active in infra-red absorption but weak 
or absent in the Raman spectrum may be the 
one involved in the resonance interaction. It is 
hoped that the infra-red spectrum may soon be 
obtained in this laboratory and that perhaps this 
additional information may make possible a more 
definite explanation of the occurrence of the 
2158 line. 
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Disubstituted Acetylenes 


The Raman data for the 2200 cm“ region of 
the spectra of several dialkylacetylenes are col- 
lected in Table III. The weak lines at 2201 and 
2280 cm~ in the spectrum of CH;C =CCH; have 
been attributed‘ to the isotopic molecule CH;C" 
=C"CH; and it seems possible that the weak 
line near 2198 cm in the spectra of the other 
dialkylacetylenes may have a similar origin 
(see IV). It seems probable that the two strong 
lines near 2233 and 2300 from Fermi 
resonance interaction between the triple-bond 
fundamental and the first overtone of the anti- 
symmetrical, infra-red active, vibration of the 
C—C=C-—C group, as for the well-known case 
of 2-butyne.” If this is true, the infra-red fre- 
quency should remain practically constant at 
1126 cm~', the value for 2-butyne. Experiments 
to test this point are a part of the future program 
of this laboratory. 

In addition to the above frequencies, a line 
appears near 2250 cm~! in the spectra of all of 
the dialkylacetylenes measured in this labora- 
tory, with the exception of 9-octadecyne and 
4-octyne. For the former, the spectrum was in- 
complete and the line may actually exist, but for 
the latter the spectrograms were excellent, with 
no trace of the line, either on long or short ex- 
posures. It is remarkable that the line should be 
absent for 4-octyne in view of its great intensity 
in the next lower member of the series, 3-hexyne. 
Several combinations of the observed Raman 
frequencies fall in this region, but it does not 
scem possible upon the basis of these combina- 


result 


TABLE II. Raman spectra of monoalkylacetylenes, 
R(CH2),C =CH, near 2200 cm=. 


Refer- 


Av(1)p ence 


2124(10) 
2118(10)0.5 2158(3)0.5 
2118(10)0.5 
2118(10)0.4 
2118(10)0.4 
2118(10)0.4" 





2068(1)P 
2064(0) 
3 2071(0) 
; 4 2064(0) 
2 2064(0) 


2094(0) 
2097(0) 
2094(0) 





* G. Glockler and H. M. Davis, J. Chem. Phys. 2, 881 (1934); Bryce 
L. Crawford, Jr., ibid. 8, 526 (1940). 

» This investigation. 

* This laboratory, see V. 
. ‘Careful remeasurement of the depolarization factor for this line 
a that the previously reported value of 0.2 should be corrected 
004, 


jisinslceepeaiiaiatmedin 


*° Cf. Bryce L. Crawford, Jr., J. Chem. Phys. 7, 555 
(1939), : 
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tion frequencies to account simultancously for 
the great intensity of the 2250 line in 3-hexyne, 
its appearance with smaller intensity in the 
spectra of other dialkylacetylenes and its strange 
absence in 4-octyne. 

The spectrum of 1-chloro-1-heptyne was ob- 
tained to see how the spectrum in the 2200 cm~! 
region would be affected by replacing the hydro- 
carbon group by a single atom of a different 

TABLE III. Raman spectra of dialkylacetylenes, 


CH;3(CH2),C =C(CH2)»,CH3, near 2200 cm™, 


m Av(])p 


2235(7)P —- 2280/1) 
2238(10) 
2238(10) 
2235(8) 
2233(10)0.4 
2233(6) 
2231(9)0.4 
2235(8)0.4 
2232(8) 
2232(8) 
2234(10)0.5 
2234(8) 
2233(8) 
2233(8)0.4 
2231(8) 
2233(8) 
2231(7)0.4 
2233(7)0.4 
2231(7) 


= 


2313(6)P 
2304(4) 
2304(4) 
2303(4) 
2305(8)0.4 
2303(4) 
2301(6)0.4 
2296(6)0.4 
mae 
2294(4 
2292(6)0.5 
2293(4) 
2292(4) 
2296(5)0.4 
2293(4) 
2294(4) 
2294(4)0.4 
2294(5)0.4 
2290(4) 





2201(1) 
2198(0) 2256(1) 


2247(8)0.4 
2249(1) 


2200(1)P 
2199(0) 


2195(0) 2245(2)0.4 


NTU WN NN HK EK KK OOOO CO | 
SUE ewer eRe Uh wre O 





8 Mean values, see Crawford's summary of data, reference 14. 


b Blanche Gredy, Théses, Paris, 1935. 


¢ This laboratory, see V. 
4 This investigation. 
e This laboratory, see ITI. 


kind. As Table I shows, the 2233, 2300 doublet 
of the alkylacetylenes is replaced by a single 
strong line at 2243 cm~. This would be under- 
standable if it were assumed that the anti- 
symmetrical, infra-red frequency near 1100 cm~! 
were shifted to a value too low to produce a 
strong resonance interaction with the triple-bond 
fundamental. The two weak lines observed at 
2148 and 2194 cm may arise from weak reso- 
nance interactions of the first overtone of the 
lowered 1100 frequency and one of the combina- 
tion frequencies that fall in this region with the 
triple-bond fundamental. 


Other Frequencies 


The strong, very broad line near 206 cm~ in 
the spectra of 1-butyne, 1-hexyne, and 3-hexyne 
unexpectedly drops to 170 cm™ in the spectra of 
the propylacetylenes 1-pentyne and 4-octyne. 
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The strong lines near 340 and 630 cm for the 
monosubstituted acetylenes and near 365 for the 
disubstituted acetylenes were observed also for 
the alkylacetylenes in V. Following Crawford," 
the frequency near 630 may be identified with 
the bending of the C=C—H angle and the ones 
near 350 may be associated with the C=C —CH; 
angle. 

All the present compounds, as well as four of 
the seven alkylacetylenes in V, have a line near 


AND R. F. 


STAMM 


500 cm~'. The present compounds also have, 
without exception, a polarized line near 800 and 
lines near 1065 and 1250 cm~. Lines near 1300 
and 1330 were found for the four compounds 
having a fairly long chain (and for five of the 
compounds in V), but not for the two containing 
only ethyl groups; these two, however, have lines 
at 1315 and 1375 cm. The doublet observed at 
1440, 1457 (mean values) for the seven alkyl- 
acetylenes in V appears here at 1436, 1457. 
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The Raman Spectra of the Sugars in the Solid State and in Solution 


I, The Raman Spectra of a- and §-d-Glucose* 


FRANK H. SPEDDING AND ROBERT F. STAMM 
Plant Chemistry Subsection, Agricultural Experiment Station} and Department of Chemistry, 
Towa State College, Ames, Iowa 


(Received December 11, 1941) 


This paper presents data on the Raman spectra of anhydrous crystalline a- and B-d-glucose 
and of equilibrium aqueous d-glucose solutions containing 5, 15, and 50 moles of water per 
mole of glucose. Ultraviolet excitation (Hg \2537A) was employed. A technique for pho- 
tographing Raman spectra of crystalline organic powders and solutions with very poor scat- 
tering power is described. It consisted essentially of the use of complementary filters of aqueous 
succinic acid (saturated, 4 mm thick) and a mercury vapor absorption cell (10 cm long, 150°C). 
Thirty-six lines are reported for a-d-glucose, 32 for the 8, and numerous new bands including 
measurements of maxima and shoulders are given for equilibrium solutions. Only the obvious 
experimental differences are discussed as the detailed discussion will be given in a later paper 


of this series. 


INTRODUCTION 


HE carbohydrates play such an important 

role in agricultural chemistry that addi- 
tional knowledge about their structures and 
chemical and physical properties is greatly to be 
desired. Because of the complexity of these com- 
pounds, accurate experimental work is difficult 
and there is a great need for the development of 
new physical methods to help the investigators in 
this field. The Iowa Agricultural Experiment 
Station, therefore, in collaboration with the 
Chemistry Department has started a program 
to study the Raman spectra of the sugars and 


* From the dissertation to be submitted by Robert F. 
Stamm to the Iowa State College in partial fulfillment of 
the requirements for the degree of Doctor of Philosophy. 

t Journal Paper No. J.954, Project 696. 


their derivatives. It is hoped that the method 
and technique of the Raman effect can be 
developed so as to be a useful tool in this field. 
If the Raman spectra of the various sugars could 
be photographed readily in the solid state and in 
solution, then from the differences in their vibra- 
tional frequencies it should be possible to identify 
the sugars by their spectra and even to determine 
quantitatively the amount of one sugar in the 
presence of another providing the mixture was 
not too complicated. Chemically this can be 
done only with great difficulty. In addition to 
this if the frequencies of a given sugar in the 
solid state differ from those in solution, some 
light might be thrown on the sugar-sugar and 
sugar-water association. Also, since the fre- 
quencies may be compared with those of mole- 
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cules thought to be similar and having better 
known geometrical structures, an insight into 
the structures of the sugars might thus be ob- 
tained. Finally, a complete set of vibrational 
frequencies would permit a dynamical analysis 
to be made with the resultant proof of geo- 
metrical structure. This last objective presents 
well-nigh insuperable difficulties if carried out 
rigorously by the usual methods. However, as 
the theory develops and the newer calculating 
machinesf{ become available, even this might be 
accomplished if the vibratory frequencies are 
available. 

In this, the initial paper of a series, we are 
reporting the Raman spectra of a- and £-d- 
glucose in the solid state and that of equilibrium 
d-glucose in solution. This aldo-hexose is one of 
the most important since it is the building block 
of both starch and cellulose. 

Several investigators of Raman spectra have 
worked with the sugars, including d-glucose, in 
the past with partial success. Whiting and 
Martin! examined a dextrose solution of 500 g 
per liter. They found no Raman lines, only an 
undiminished continuum. Kutzner? obtained 
three lines of very low intensity accompanied by 
a heavy continuum after 30 hours’ exposure. 
Wiemann’ used charcoal to decolorize a 30 per- 
cent dextrose solution which he then filtered 
through sintered glass and reported that very 
little continuum was present in his pictures. He 
found eleven bands in solution. A. L. Sundara 
Rao‘ reported twelve lines for d-glucose solution, 
two of them broad. Only the last investigator 
reported any lines in the C—H and O—H regions 
and he stated that the spectra were so faint that 
it was difficult to distinguish them from the 
background. 

Since the glucose frequencies can also be ob- 
served in infra-red absorption, it should be 
mentioned that Coblentz® studied molten dex- 


} Those used for evaluating secular determinants. 

' Whiting and Martin, Trans. Roy. Soc. Can. III [3] 25, 
Sites) (1931). (Visible Hg excitation, quinine sulphate 

ter. 

*Kutzner, Naturwiss. 20, 331-332 (1932). (Super- 
= solution at 20°C, quartz Hg arc, visible excita- 
ion. 

*Wiemann, Comptes rendus 203, 789-791 (1936). 
(Visible Hg excitation.) 

‘A. L. Sundara Rao, Indian J. Phys. 14, 207 (1940). 
(Visible Hg excitation.) 

* W. Coblentz, Carnegie Inst. of Wash. Pub. 65 (1906). 


trose and obtained a band from 3—-4y and a 
maximum at 4.8u. Rogers and Williams® used 
saturated solutions 0.002 mm thick and inves- 
tigated the region 2-124. They reported seven 
bands but could not resolve details in the C—H 
and O—H regions. Barr and Chrisman’ evapo- 
rated saturated solutions on cover glasses and 
used the resulting syrups. Two bands were found 
by them for each of a series of sugars including 
d-glucose. 

The phenomenon of muta-rotation leads to the 
conclusion that in aqueous glucose solutions both 
the a- and §-forms of the pyranose ring are 
present in equilibrium with one another and 
perhaps in equilibrium with other forms. There- 
fore any spectrum of the solution will be com- 
posed of a superposition of the spectra of all the 
various forms present and many of the bands 
observed will in reality be envelopes of several 
lines or bands from these forms. Since all of the 
above work (that of Coblentz excepted) was 
done with equilibrium solutions, the values of 
the vibrational frequencies obtained by them 
should be compared with our solution values in 
Table I. 

The a- and 8-forms of d-glucose are geometri- 
cally different and should possess different vibra- 
tional frequencies. It is therefore desirable to 
study the spectra of the pure forms alone. This 
can be done with the sugars in the solid state. 


EXPERIMENTAL 


From the workof Kwiencinskiand Marchlewski*® 
and especially from that of V. Henri and 
Schou? it was decided that pure d-glucose does not 
absorb light of wave-lengths longer than 2300A 
and that the 2537 mercury line could therefore 
be used as a primary exciter of the Raman lines. 
We employed one of the low pressure argon- 
mercury arcs.'® Unfortunately the direct light 
from this arc decomposes all sugars so that all 
light below 2300A has to be removed by a filter. 


6 Rogers and Williams, J. Am. Chem. Soc. 60, 2619-2621 
(1938). 

7 Barr and Chrisman, J. Chem. Phys. 8, 51-55 (1940). 

8 Kwiencinski and Marchlewski, Bull. Soc. Chim. 43, 
725-743 (1928). 

9V. Henri and Schou, Zeits. physiol. Chem. 174, 295 
(1928). 

10 The Hanovia Chemical and Manufacturing Company, 
Newark, New Jersey (Helicoidal type, 0.1 amp. at 4400 v). 
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5 moles 
H2O: mole 


Anhydrous 
glucose 


B-d-glucose Int. 


Anhydrous 
a-d-glucose Int. 


142 1 137 1 96) 
168 1 163! 1 i. 171] 
224 1 247) 
288 
poo 
316) 
322) 
345 | 
368) 


386) 
445 | 
467 
483) 
502 | 
518! 
540| 
550} 
562 | 


584| 


589} 
628 
640 
652 


702) 
715} 


1020 
1038! 


1057 
1072 


1109 


This was done by circulating a saturated (30°C) 
aqueous solution of succinic acid continuously 
through a fused-quartz filter jacket interposed 
between the arc and the Raman tube. The 
thickness of the filter solution was 4 mm and it 
was circulated by means of a midget centrifugal 
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HO: mole 
glucose 
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50 moles 
H.O: mole 
glucose Int. 


Previous investigators 
solutions only 


21) 90) 
87 | i. e72 
52) 255) 
255 
290| 
325 
336) 
349 | 
363 


327) 


389) 
404 | 
421 | 

. 447| 

? 464) 
471} 
486 | 
497 d 
520) a 519 +11(W), 513(K) 
541} : 5493(R) 


423(R) 
435 £25(W), 454(8) 


559| 
573) 
586) 
609 | 


6241(R) 


8$56(R), 853 £10(W) 


919(R), 911 £14(W) 


1000(RW) 
1020 +10(W) 


1064; br. (R), 1067 +12(W) 


1100(RW) 


115910 br. (R), 1145+9(W), 
1140(K) 


1200(W) 


1270(R), 1262 +8(W) 


pump." A number of other solutions were tried 
but proved unsatisfactory ; ¢.g. sodium succinate 
tended to give off gas which permitted the short 
ultraviolet light to get through where the bubbles 


1! Obtainable from the Eastern Manufacturing Company, 
New Haven, Connecticut. 
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TABLE I—Continued. 


5 moles 
HO: mole 
glucose 


Anhydrous 
B-d-glucose 


Anhydrous 
a-d-glucose Int. 


1295 
m. 1331 5 


Int. Int. 


m. 1370} 


m, 


1507 
1605 | 

i. 1646) 
1688 
2057 

i. 2168 
2280 
2447 
2528 
2549) 
2556) 


vg(H20) 


(HO) 


br.m. 2625 


2640) 


br. {2706 
i.flat{ 2744 
m. {2781 | 


sh. 2801 


' 
br.m. 2896) 
2916 } 


br.m, 2945| 
3022} 


2946 
2961 


| 3237 .| 


(3349) 
br.m.4 3439 
3522| 


22| 


3396 


3503} 
3556} 
3714| 


15 moles 
1120: mole 
glucose 


10 m. 
sh. 3 m. 


. 1460 


i. 2140) i. 


. 2902 


. 2943 


. 3434! 
* 3619] 


50 moles 
HO: mole 


Previous investigators 
glucose 


solutions only 
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| Observed on only one plate. 


2 From the relative intensities of Hg \2563.9 and 2578.4A, we conclude that this maximum is not the mercury line 2563.9A but is a Raman line. 
br. broad, c. center, d. diffuse, i. center interpolated, m. maximum, sh. shoulder, ? doubtful. (K) Kutzner, reference 2, (R) Rao, reference 4 (sub 
scripts are his intensities), (RW) Rogers and Williams, reference 6, (W) Wiemann, reference 3. 


formed and the sugar was then decomposed 
because the bubbles did not always leave the cell 
immediately. The succinic acid’ proved very 
satisfactory. The ultraviolet arc emits the 1850A 
line of mercury with considerable intensity. This 
in turn generates ozone in the air which tends to 
cut down the intensity of ultraviolet light reach- 
ing the spectrograph. The difficulty was over- 
come by having a strong suction just above the 
arc so that the ozone was removed as soon as it 
was formed. 

The sugar crystals were contained in a fused- 
quartz tube of 6 mm inside diameter with a 
fused-quartz window sealed on the end. The 
length of the sample was about 5 mm. Longer 


®" The quantitative absorption data can be found in the 
oper by Ley and Arends, Zeits. f. physik. Chemie B17, 
77 (1932). 


samples were not needed since very little light 
penetrated to the window from distances greater 
than 5 mm. In order to prevent light which had 
not penetrated the crystals from reaching the 
spectrograph a brass collimator containing two 
diaphragms of 4-mm aperture placed 3 inches 
apart was affixed to the tube with Duco cement. 
The interior of the collimator was painted a flat 
black. To gain better light transmission some of 
the very fine crystalline powders were moistened 
with absolute methyl alcohol. Only one line of 
the methyl alcohol appeared on any of the plates 
(Av 2836 cm). With the dispersion used, no 
change in the frequency of the sugar vibrations 
could be detected when this procedure was used 
and the exposure time was cut frequently by a 
third. 
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The sugar solutions were filtered through 
sintered glass and placed in a fused quartz 
Wood's tube with a flat, polished window." This 
tube had a diameter of 15 mm and a possible 
active length of 5 inches. A half-inch length of 
the front of the tube and the horn were painted 
a flat black thus leaving a 3-inch effective length 
for irradiation. The solution tube had a col- 
limating tube also and was used in conjunction 
with the filter jacket mentioned above. The 
mercury arc was surrounded with a polished 
sheet aluminum reflector, and the entire as- 
sembly was located in a large box. The collimator 
tube projected through the box, thus eliminating 
the possibility of light from the arc getting into 
the spectrograph. In the case:of powders there 
can be no (or at most only a very short) dead 
optical space in the front end of the sample, for 
unless the powder near the window is illuminated, 
the secondary radiation cannot get through to 
be photographed. Thus with the powders con- 
siderable primary light (except 2537 which was 
especially removed) got into the spectrograph 
due to multiple reflections from the crystal faces, 
and the continuum on each side of 2537 was 
strong as were the rest of the mercy lines. Since 
this arc yields a particularly strong resonance 
line (the makers claim it to be nearly 80 percent 
monochromatic 2537), on long exposures, we 
found that the usual practice of having mercury 
vapor in the spectrograph to absorb the 2537 
line was not efficient enough and we therefore 
employed a mercury vapor absorption cell. This 
cell was electrically heated, had an absorbing 
column 10 cm long, and when heated to ca. 150°C 
with a drop of mercury in it, completely ob- 
literated the 2537 line. 

The spectrograph employed was a Bausch and 
Lomb medium quartz type. In the region used 
the linear dispersion ran from around 73—-9A/mm. 
The room housing the spectrograph was ther- 
mostated to within 1°C. 

The sugars used were anhydrous a-d-glucose," 
[a]p**°=110.4 initial rotation (Lit. 111.5); 
[a]p**°=52.2 final rotation (Lit. 52.7): a-d- 


13 This and the filter jacket were obtained from The 
American Instrument Company. 

4 “Cerelose” from the Corn Products Company, New 
York. 
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glucose monohydrate,'® [a ]p"?©=100.8 initial 
rotation (Lit. 101.4); [a ]p*°©=47.5 final rota- 
tion (Lit. 47.9): anhydrous B-d-glucose (two 
samples), one (from The Pfanstiehl Chemical 
Company), [a ]p**=18.5 initial rotation (Lit. 
18.4); [a]p*°©=52.5 final rotation: another 
(prepared according to Hudson and Dale'®), 
[a]p**=18.5 initial rotation; [a]p*°=51.8 
final rotation. The literature specific rotation 
values given were taken from the paper by Isbell 
and Pigman.!” Unfortunately their figures do not 
agree exactly with those of other investigators 
and the discrepancies are of the same magnitude 
as those between our values and those of Isbell 
and Pigman. Therefore, all we can say about the 
purity of our a- and 8-sugars is thatany impurities 
or other forms present must be there in such 
small amounts that their Raman spectra would 
not be observed. The solution pictures were made 
with d-glucose purified by the method of Hudson 
and Dale’ using ‘‘Cerelose”’ as starting material." 
A solution of ‘‘Cerelose’’ is usable but is slightly 
yellow originally and is a very bright yellow 
after a 24-hr. exposure. Also the impurities which 
can be removed by clarification and recrystal- 
lization increase the exposure time and give rise 
to fluorescence thus tending to hide the fainter 
Raman lines. 

It was found difficult to measure the lines 
under our comparator since due to its magni- 
fication and to the low contrast of the fast plates 
used, the lines tended to wash out. Therefore the 
spectrograms were enlarged about sevenfold on 
high contrast paper and the wave-lengths of the 
Raman lines found by the Hartmann formula. 
Such was not the case, however, for the solution 
bands since they were frequently envelopes of a 
number of lines often diffuse and ill defined. It 
was found better to microphotometer these 
plates with a tenfold enlargement and to measure 
the band edges, shoulders, and maxima on the 
curves and to determine their wave-lengths by 
means of the Hartmann formula. Checks of the 
solution measurements were made on enlarge- 


6 Prepared according to Hudson and Yanovsky, J. Am. 
Chem. Soc. 39, 1013 (1917). 

16 Hudson and Dale, J. Am. Chem. Soc. 39, 320 (1917). 

17 Isbell and Pigman, J. Research Nat. Bur. Stand. 18, 
158 (1937). 

18 In the clarification both ‘“‘Norit’”’ and SOz were used 
prior to the concentration and subsequent crystallization. 
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Fic. 1. Enlargements representative of the Raman spectra obtained by Hg \2537A excitation. 
1. Anhydrous a-d-glucose. 2. Anhydrous B-d-glucose. 3. 5 moles H2O/mole glucose. 4. 15 moles 
H20/mole glucose. 5. ca. 520 moles H2O/mole glucose (0.1 molar). The two bands in 1 and the higher 
of these two bands in 2 lying to the low side of Hg 2803.5A belong to a band system in the arc 


which encroaches slightly on the O-H region. 


ments to be sure that no false information (due to 
scratches or flaws on the plates) would be con- 
veyed by the curves. Too much stress must not 
be put on these band edges since their position 
varied slightly with the time of exposure. The 
low frequency ‘‘crystal” lines on the crystal 
plates were measured on the comparator after 
first pricking the lines with a fine needle.'® In 
general the measurements are judged to be 
accurate to within +2 cm~'. The intensities are 
qualitative only and are based on a scale of 0 to 
20 with 2893, 2895-2908, and the 2900-cm— band 
of the a-, B-, and solution plates respectively 
being given an arbitrary rating of 20. Plates were 
made of the direct arc in order to be sure which 
lines originated in the arc and which in the sugar. 

Below is given a resumé of plates measured. 
After the number of the plate, the words “‘dry”’ 
or “‘moistened”’ indicate whether the powder was 
photographed dry or moistened with absolute 
methyl alcohol. Next is given the slit width on 





Some of these low Ay lines seen in the reproductions 
may arise in the arc. They are Av 47, 79, 111, 191, and 380 
cm™. These arc lines are sharp but faint and probably too 
many crystal frequencies were rejected by being assigned 
to these weak lines in the source. We have leaned over 
backwards in this respect. 





the spectrograph and finally the number of 
hours which the plate®® was exposed. 

a-d-glucose anhydrous (all I-O plates): (1) 
dry, 0.025 mm, 35 hr. (2) dry, 0.02, 25 (3) 
moistened, 0.02, 15 (4) moistened, 0.02, 25. 

B-d-glucose anhydrous (all I-O plates): (1) 
moistened, 0.02, 24 (2) dry, 0.02, 60 (3) dry, 0.04, 
72 (4) moistened, 0.03, 40 (5) dry, 0.03, 42 (6) 
dry, 0.03, 72. 

Saturated equilibrium d-glucose (ca. 5 moles 
H,O/mole glucose) : (1) I-O, 0.025, 25 (2) 103-0, 
0.04, 75 (lines also confirmed on two other plates 
not measured). 

Dilute equilibrium d-glucose (ca. 15 moles 
H,O/mole glucose) : (1) 103—-O, 0.03, 84. 

Dilute equilibrium d-glucose (ca. 50 moles 
H,O/mole glucose): (1) 103-O, 0.03, 47 (lines 
confirmed on other plates not measured). 

The results are summarized in Table I. The 
strong lines were observed and measured on all 
plates mentioned. In the first three columns the 
weaker lines were observed and measured on at 
least two plates unless marked by a superscript 
1 in the table. 


20 All Eastman plates. 
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Fic. 2. Comparison of the Raman frequencies (in cm~) of anhydrous a- and 8-d-glucose. Intensity estimates 
plotted vertically. 


DISCUSSION OF RESULTS 


In the crystals there is only a fair one-to-one 
correspondence among the vibrational frequen- 
cies of a- and 6-d-glucose below 1500 cm~. One 
notable discrepancy is the absence of a strong 
frequency in the @-d-glucose near the rather 
intense 842-cm™ line of the a-form. While it is 
a fact that a strong frequency exists in this 
region for a whole series of six membered hetero- 
cyclic ring compounds and corresponds to the 
breathing frequency, we do not wish to make any 
conjecture about the absence of a strong fre- 
quency of the 8-form in this region until more 
intense spectrograms have been taken along with 
accurate intensity measurements and _ polariza- 
tion characteristics of the lines. We hope to carry 
out these experiments if we are successful in 
growing large single crystals free from fluorescent 


impurities and flaws. The most striking difference 
between the two crystalline forms is exhibited in 
the 2900-cm- region where the form of sugar 
may be identified by a glance at the spectra. The 
O—H frequency region is nearly bare for both 
forms with only one line in the a and one band 
in the 8. On the best plates, which we have 
reproduced in Fig. 1, there is some evidence in 
the @ for another frequency of slightly smaller 
value than the O—H frequency mentioned. 
Rough preliminary measurements on the a-mono- 
hydrate seem to show no marked differences 
from the anhydrous a. This seems astounding 
and we are investigating the spectra further. We 
do not feel that speculation as to the identity of 
the various frequencies or as to their appearance 
or non-appearance would be very profitable at 
this time. We are, therefore, reserving any such 
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Fic. 3. Superposition of microphotometer tracings of: 1. 5 moles H»O/mole glucose, 2. 15 moles HzO/mole glucose, 
3. 50 moles H,O/mole glucose. 4. High part of plate with suitable exposure in this region but under-exposed in lower \ 


region. 5 moles H2O/mole glucose. 
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comments until we have obtained much more 
data on the various sugars. A comparison of the 
Raman frequencies of anhydrous a- and £-d- 
glucose is given in Fig. 2. 


Glucose Solutions 


When the sugars are dissolved in water, 
allowed to come to equilibrium, and the Raman 
spectra of the solutions obtained, the frequencies 
are no longer sharp but are broadened consider- 
ably. Below 1500 cm~ there are 15 solution 
bands having no counterpart in the a- or B-form. 
There are no a-lines which are not covered by a 
solution band and only two 6-lines (Av 825 and 
1507 cm~'). However, there are four a-lines 
(Av 363, 887, 1057, and 1345) lying more than 
5 cm from band maxima or shoulders in all 
three concentrations and five B-lines (Av 1051, 
1118, 1312, 1363, and 1452) in this category. 
Unfortunately, the microphotometer available to 
us was not very satisfactory in that there was a 
certain amount of random drift in the intensity 
deflection. While this could be kept at a minimum 
in closely adjacent regions, absolute intensities 
or true shapes of lines were not to be depended 
on except as rough estimates. It is hoped that 
this microphotometer will be in first class shape 
in the future as soon as parts can be obtained. 
In Fig. 3 we have reproduced the microphotom- 
eter curves of the three solutions of different 
concentrations. Within the accuracy of our 
tracings, which we estimate to be of the order of 
magnitude of 5 cm~', there are no systematic 
changes with dilution except for the frequency 
894-4899-—905. In order to obtain reasonable 
intensities for the frequencies below 1500 cm! 
it is necessary to over-expose’ the 2800-3600 cm7! 
region. This has a tendency to flatten the tracing 
in that range. In curve 4 we have micropho- 
tometered a plate (5 moles H2,O/mole glucose) 
where the 2800-3600 cm~! region was not over- 
exposed. 

We do not yet feel justified in drawing any 
conclusions from the curves as to the percentages 
of a- and 6- (or other) forms present in solution. 
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A careful study of how the frequencies and inten- 
sities of the lines shift in going from the an- 
hydrous crystalline forms to the various crystal- 
line hydrates (if they exist) should indicate how 
the frequencies are affected as the sugar-sugar 
association is replaced by the sugar-water 
association. More information should also be 
obtained by studying the changes in the shape 
and position of the line envelopes in the equi- 
librium sugar solutions as the concentration 
changes from a saturated solution (where con- 
siderable sugar-sugar association must still exist) 
to a very dilute solution (where only sugar-water 
association could be present). From the above 
work it can be seen that these shifts must be 
small, about the order of magnitude of our errors. 
The question as to whether the difference in the 
sugar-sugar sugar-water association can cause 
enough changes in the frequencies and intensities 
to permit the solution bands to be completely 
accounted for by the a- and 8-forms will have to 
wait until we have photographed more sugars 
and their hydrates and have obtained precise 
microphotometer tracings of their solutions. 


The Water Bands 


These are located*! at ~ 200, 320-1020, 1650, 
2167, 3440, and ~ 4000 cm~. All of them but the 
last may be seen in the most dilute d-glucose 
solution (0.1 molar), reproduced in the fifth 
picture in Fig. 1. Although the water bands in 
this case cover almost every region except 1000- 
1500 cm it can be seen from this last picture 
that even in dilute solutions the Raman lines of 
the sugar can be seen superimposed upon the 
water bands. 
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The Henry’s law constants for HT evaporating from approximately 10~°M solutions in 
normal He and for DT from nD» were determined by distillation and measuring the radio- 
activities of the distillate gases. On the basis of Lewis’ and Hanson’s data showing Hz and D2 


to form nearly perfect solutions, the Henry’s law constants were taken to represent the vapor 
pressures of the pure HT and DT liquids. A plot of vapor pressures of various kinds of hydrogen 


at 20.4°K led to 45+10 mm of Hg as the predicted vapor pressure of liquid nT, at 20.4°K. 





HE vapor pressures of the hydrogen mole- 

cules with tritium (radioactive hydrogen)! 
atoms should help settle certain points in the 
theory of liquid hydrogen involving the mass and 
zero-point properties of the liquid. In addition, 
a distillation technique may prove to be desirable 
for the concentration of tritium. 

The method used consisted principally of the 
distillation under equilibrium conditions of a 
small volume of the liquid containing some of 
the tritium-containing molecules, samples of the 
gas being taken frequently and their specific 
radioactivities determined. Under such condi- 
tions of no refluxing and perfect stirring of the 
liquid, the Rayleigh? distillation formula applies, 
with a slight modification due to the large frac- 
tion of the hydrogen in the gas phase. 

Suppose the vapor pressure of HT over the 
very dilute solutions (approximately 10-°M) in- 
volved here is given by 


Pur=B6Nur, 
where Pyr is the pressure and Nur is the mole 
fraction of HT, and, similarly, 
Py,=aNy,=a(1— Nur). 
Then the specific radioactivity or mole fraction 
of HT of the distillate gas p will be given by 
cee ’) (“Da 
= po 
1+ (a/B)v/V 


where V is the volume of the gas space over the 
liquid and v» and v are the volumes of the liquid 
at the beginning of the distillation and at any 

'Luis W. Alvarez and Robert Cornog, Phys. Rev. 56, 


613 (1939). 
2 Lord Rayleigh, Phil. Mag. [5] 4, 527 (1902). 





subsequent point, measured in terms of the 
volumes of gas to which they are equivalent. 

To derive this expression one notes that the 
total number of moles of HT in the distillation 
vessel is 

t= pV+ p(a/B)u, 

since p is the mole fraction of HT in the gas and 
is equal to (8/a)Nur. 

For a small removal of gas from the distillation 
chamber, changing v by dz, 


dt= Vdp+(a/8) pdv+(a/B)vdp 
or, after dividing by dv and replacing dt/dv by p, 
p=(V+(a/B)v)dp/dv+a/Bp. 
Rearranging gives 
(8/a—1)d In (V+(a/8)v) =d In p, 


since V is essentially constant. Finally, integra- 
tion gives the equation desired. Figure 1 illus- 
trates the expected course of the specific activity 
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Fic. 1. Activity distillation curve for HT from mH: 
(20.4°K). Solid line: calculated curve for pressure ratio o! 
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295°K and one atmosphere) and a/8 was 3.0. 
The distillation bulb was immersed deeply in 
a bath of boiling normal H: and the previously 
prepared radioactive gas introduced with a 
Toepler pump through several feet of small 
bore (1 mm) capillary glass and copper tubing. 
After a few minutes, gas was slowly withdrawn 
at equilibrium pressure into a series of bulbs 
which were stored for later measurement of their 
radioactivities. The total volume of gas used 
varied between one and five liters NTP and the 
sample bulbs were mainly of 25 cc and 100 cc 
capacities, the smaller bulbs being used near the 
point of disappearance of the last of the liquid. 
Obviously under the conditions outlined reflux- 
ing of the distillate is quite impossible. The 
second condition, that the liquid be in equi- 
librium with the vapor, is less easily established. 
Indeed, for the run shown in Fig. 1 some five 
liters of gas were used and only the last liter 
distilled under equilibrium conditions as attested 
by the agreement between the points and the 
calculated curve. The result is understandable 
when one realizes that the vessel in which the gas 
and liquid came to equilibrium was about 2 


TABLE I. 








Vapor pressure at 20.4°K 





Molecule (mm Hg) Source 
H, 760 
HD 438 Reference 1. 
D, 256 References 2-5. 
HT 254+16 This research. 
DT 123+6 This research. 
To (45+10) Extrapolation, Fig. 3. 








'R. B. Scott and F. G. Brickwedde, Phys. Rev. 48, 483 (1935). 

2R. B. Scott, F. G. Brickwedde, H. C. Urey, and M. H. Wahl, 
J. Phys. Chem. 2, 454 (1934). 
sas). Clusius and E. Bartholomé, Zeits. f. physik. Chemie B30, 237 


‘E. Bartholomé, Zeits. f. physik. Chemie B33, 387 (1936). 
*H. D. Megaw and F. Simon, Nature 138, 244 (1936). 


square cm? in cross-sectional area and about 15 
cm long. Five cc of liquid hydrogen created quite 
a distance for diffusion to accomplish stirring. 
Fortunately, diffusion was sufficiently rapid when 
the liquid volume was reduced. 

Figure 2 presents similar data for distillation 
of DT out of nD». In this case the bulb size was 
somewhat greater, reducing the number of points. 
However, the larger size of the bulbs increased 
the accuracy of the counting because the 500 cc 


VAPOR PRESSURE 


for the case where V was 500 cc (measured at 
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Fic. 2. Activity distillation curve for DT from nD,» 
pay Solid line: calculated curve for pressure ratio of 
1 to 2.08. 


counter used could be filled several times to the 
2 cm of Hg pressure used. The point placed just 
at the point of disappearance of the liquid was 
obtained by using the activity of the gas present 
after the liquid had disappeared. The gas volume 
V was 550 cc in this run. 

The HT samples were prepared by allowing Li 
metal which had been irradiated for several 
months by the Berkeley cyclotron (the authors 
are indebted to the Radiation Laboratory and 
Dr. Robert Cornog for this irradiation) to react 
with H,O after, all air had been removed by 
pumping on the liquid-air-cooled reaction mix- 
ture. The appropriate amount of Li metal was 
cut from the irradiated block and placed in a 
liquid-air-cooled tube containing sufficient water 
as ice. The tube was then placed on the vacuum 
line and pumped. The vessel then warmed to 
room temperature and evolved He. The gas was 
passed through liquid-air traps for water re- 
moval. The procedure for DT was identical, 
except that D.O was used instead of HO for 
reaction with the lithium. The gases were stored 
in large bulbs and transferred with a Toepler 
pump to the distillation cells. The purity of the 
gases must have been high because the capillary 
tubing immersed in liquid H2 for about two feet 
would have plugged immediately with any air, 
water vapor or any gas except the lightest rare 
gases. In addition the Hz and D2 samples gave 
the vapor pressures reported in the literature. 

Two entirely separate runs were made with 
HT. The first one gave a value for a/8 of 3.34.5 
considerably less accurate than the one de- 
scribed above, but in definite agreement with it. 
There seem to be only two likely sources of 
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Fic. 3. Plot of vapor pressures of hydrogen molecules 
versus mass, at 20.4°K. 


error in this type of experiment. The first is the 
condition that diffusion equilibrium exist and 
the second is that the samples of gas have their 
radioactivities measured properly. In connection 
with the first point, one notices on Fig. 1 that 
there remained only 600 cc of hydrogen as liquid 
when the data were considered to represent safely 
equilibrium distillation. When one realizes that 
this amounts to about 3 cc of liquid hydrogen 
in the distillation cell or a layer of liquid hydro- 
gen approximately } cc thick over the 2 cm? 
bottom of the cell, it is difficult to imagine any 
but equilibrium conditions applying after 10 
minutes, which was the average interval between 
samples. 

The gaseous samples were counted by filling 
the counter with approximately 1.5 cm of Hg of 
He or Dz and 2 or 3 mm of Hg of ethyl alcohol 
vapor. A standard uranium sample was used to 
check the reproducibility of the counting appa- 
ratus. In connection with the accuracy of the 
counting, it was observed that air had not 
entered the sample bulbs and that the gas com- 
position was uniform by noting the operating 
voltage for the counter when filled with the 
various samples. Counters are very sensitive in 
this respect. This being so, one could feel certain 
that the counter’s sensitivity to tritium radiation 
would be reproduced if its sensitivity to uranium 
were reproduced. 

The figures quoted for the errors in the values 
of the vapor pressures were obtained as shown 
in Fig. 1 by sketching in the curves for values 
of a/8 such that most of the points were included 





between the limits. The values for V were 
obtained by measurement of the cell volume and 
correcting to 20.4°K. 

Table I summarizes the results for the vapor 
pressures of the hydrogen molecules at 20.4°K. 

In compiling the data the Henry’s law con- 
stants 8 obtained for the very dilute solutions 
were assumed to represent the true vapor pres- 
sure of the pure solute, i.e., the solutions were 
assumed to be perfect. Evidence for the essential 
correctness of this assumption exists in the data 
on the vapor pressures of He, De solutions.’ 

The first point the data emphasize is that the 
vapor pressure of HT is essentially the same as 
that for normal Ds». This is rather to be expected 
when one realizes that the difference in moment 
of inertia can have little effect so long as nearly 
free rotation occurs in the liquid, a condition 
fairly well substantiated by specific heat data 
on Hao. 

The second point is shown in Fig. 3 where the 
vapor pressures of the various hydrogen mole- 
cules are plotted against their masses, the value 
for T. being obtained by extrapolation of the 
curve. Very evidently, the distillation technique 
should be successful in concentrating tritium if 
one were to keep the equilibria 


2HT=T.+H, 
and 
2DT=T.+D., 


well catalyzed by periodically passing the gas 
over a heated catalyst in the last stages of the 
concentration. 

The complete significance of the curve in Fig. 3 
must await a successful theory of liquid hy- 
drogen. However, qualitatively it reveals the 
effect of translational! zero-point energy on the 
free energy of the liquid. The van der Waals 
forces holding the liquids together will be nearly 
identical for all the liquid hydrogens. Therefore 
one must expect the mass effect observed to be 
associated very directly with the differences of 
the zero-point energies of the compressional 
vibrations of the several liquids. The effects due 
to rotation certainly are small compared with 
these mass effects. In addition they are averaged 
out somewhat in the data chosen for Fig. 3. 


3 Gilbert N. Lewis and W. T. Hanson, Jr., J. Am. Chem. 
Soc. 56, 1000 (1934). 
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X-ray studies of solid solutions of alkali halide pairs with a common ion at room temperature 


and at 550°C show that these pairs can be divided into three groups depending on the percent 
deviation 6 of their lattice parameters. For values of 6 less than 6 percent miscibility is com- 
plete at room temperature. For values of 6 between 6 and 13 percent miscibility is complete at 
550°C. For values of 6 larger than about 13 percent miscibility is not complete at 550°C. In 
terms of the Born theory of ionic lattices, the free energy of mixing, and Vegard’s law, an 
approximate theory for the phase curve of these solid solutions is derived. For a typical alkali 
halide with an electrostatic energy of 180 kcal. at the absolute zero of temperature the equation 


of the phase curve is 


98/T(1—2x)+I1n x—In (1—x) =0, 


where x is the mole fraction of one component dissolved in the other. The temperature above 
which miscibility is complete is given by T=4.56. The heat of mixing is approximately 
—98* cal. The predictions of this theory are in good agreement with the above-mentioned 


experimental observations. 


SOLID SOLUTIONS OF THE ALKALI HALIDES 


T has long been known that when certain 
mixtures of alkali halides with a common ion 
are brought into contact solid solution formation 
takes place.! Only recently, however, has the rate 
of attainment of equilibrium for this solid solu- 
tion formation been carefully studied.? It has 
been found that at room temperature if the salt 
mixtures are kept dry, it will take a matter of 
several weeks for equilibrium to be reached. The 
presence of moisture will cut thé time to one day. 
At high temperatures if the salt mixtures are 
fused, quickly quenched, then ground in an agate 
mortar and annealed in a Victor Meyer furnace 
at 550°C, it will take only a few hours for equi- 
librium to be reached. 

Continuing the work of Matsen and Beach in 
this laboratory, a study was made to determine 
what are the stable phases when equimolal 
mixtures of alkali halide pairs with a common 
ion are brought together at room temperature 
and at 550°C. In general, each component is 
somewhat soluble in the other and two solid 
solution phases occur. The extent of solution is 
greater if the lattice parameters are more nearly 
equal, and the amount of solubility increases 
with the temperature. For certain mixtures for 





wong a general review see Bruni, Chem. Rev. 1, 345-375 


. 


oe and Beach, J. Am. Chem. Soc. 63, 3470-3473 
ss 
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which the lattice parameters are nearly equal 
miscibility is complete at room temperature and 
a single phase will result when equilibrium is 
reached. For others, for which the percent dif- 
ference between the lattice constants is somewhat 
larger two phases result at room temperature but 
miscibility is complete at 550°C. For a third class 
of these mixtures for which the lattice parameters 
are widely different the solubility will be small 
and two phases will result for all temperatures 
below the melting temperature. 

The method used in this study to determine 
the equilibrium phases is an x-ray powder 
method fully described in the paper by Matsen 
and Beach. No attempt was made to determine 
quantitatively the composition of the phases, 
but merely to determine whether or not mis- 
cibility was complete. In most cases where it 
was found that two phases were formed, these 
phases were found to be nearly pure phases on 
the criterion that their lattice constants were not 
much different from the lattice parameters of the 
pure phases. A similar study had previously been 
made by the x-ray method for some of these 
mixtures at room temperature.* Attempts to 
determine the width of the miscibility gap of the 
phase diagram had also been made by other 
methods.' The various results have been collected 


3 Havighurst, Mack, and Blake, J. Am. Chem. Soc. 47, 
pp. 29-43. 

















































PI 







eer als 












Alkali halide 
pairs 





ARTHUR V. 
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TABLE I. 
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Remarks from other studies 








NaF-KF 
NaClI-KCl 


NaBr-KBr 
Nal-KI 


NaF-Nal 
NaCl-NaBr 
NaBr-Nal 


KF-KCl 
KF-KBr 
KF-KI 


KCI-KBr 
KCI-KI 


KBr-kI 


10.8 
9.1 


39.8 


8.8 


17.8 
23.5 
32.3 


4.8 
12.3 


7.1 


15°C Two nearly pure phases* 
550°C Two nearly pure phases 

15°C Two nearly pure phases 
550°C Complete miscibility 


15°C Two phases 

550°C Complete miscibility 
15°C Two phases 

550°C Completely miscible 


15°C Two nearly pure phases 
550°C Two nearly pure phases 
15°C Complete miscibility 
550°C Complete miscibility 
15°C The x-ray lines could not 
be interpreted 
550°C Complete miscibility 
15°C Two nearly pure phases 
550°C Two nearly pure phases 
15°C Two nearly pure phases 
550°C Two nearly pure phases 
550°C Believed to be two nearly 
pure phases. KF lines do 
not show very well 
15°C Completely miscible 
550°C Completely miscible 
15°C Slight dissolving of KCI 


in 
550°C More extensive dissolving 
ing of KCI in KI but 
not completely miscible 
15°C Two phases 
550°C Completely miscible 


H.M.B.:+ Solution of a few percent KCl in NaCl 

Bruni:{ Mixed crystals formed from melt destroyed 
at room temperature 

Bruni: Mixed crystals formed from melt destroyed 
at room temperature 

H.M.B.: Very low solubility at room temperature 

Bruni: Mixed crystals formed from melt destroyed 
at room temperature 


Studied by many investigators 


Bruni: At room temperature the miscibility gap 
extends from 7-90%. Near the melting tempera- 
ture the miscibility gap is between 49 and 91% 


H.M.B.: KBr soluble to the extent of 15 or 20% in 
KI at room temperature 
Bruni: Miscibility gap between 34 and 93% 








* By this phrase we mean that the x-ray picture shows two phases that differ very slightly from the pure components. 
t Reference 3, Havighurst, Mack, and Blake. 


¢ Reference 1. Within the approximations of the theory to be developed 100(a2—<a1) /a: = 100(a2 —a1) /a2 =100( a2 —a1)/ [(a1 +2) /2] =5, 
that it does not matter which particular definition we adopt for purposes of tabulation. 


in Table I. The quantity 6 that appears in this 


table is the percent deviation between the lattice 
constants a; and a2 of the component salts. 
It will be seen, as had been noted previously,* 





that for salt mixtures for which the value of 6 
is less than 5 or 6 percent miscibility is complete 
at room temperature. Working under the super- 
vision of Dr. Beach in this laboratory I have 
made further x-ray studies at high temperatures 
and noted that for values of 6 between approx- 
imately 6 percent and 13 percent miscibility is 
partial at room temperature and complete at 
550°C (i.e. only one phase is observed in the 
x-ray powder photograph). For larger values of 6 
two phases are always observed and especially 
at room temperature the phases are very nearly 
the original pure components. 

In order to explain these phenomena, I have 
developed a theory for the phase diagram of 


<0 





these salt mixtures in terms of the Born theory 
for ionic lattices, Vegard’s law, and the free 
energy of mixing. It is commonly accepted that 
in a solid solution such as KCI-KBr the Cl ions 
substitute for Br~ ions at random in the KBr 
lattice. Superstructure lines which characterize 
ordered lattices such as occur in many alloys 
have never been observed for mixed crystals of 
alkali halides. One must therefore assume that 
if there is any ordering phenomena the Curic 
point for order-disorder occurs far below room 
temperature. In this theory it has been assumed 
that the substitution is perfectly random at all 
temperatures, i.e. whatever mixed crystals occur 
are perfect mixed crystals. 

A general expression for the lattice parameter of 
a solid solution in terms of the parameters of the 
components has been suggested to be of the form 


a" =a;"+ (a2" —a4")x, 
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SOLID SOLUTIONS OF 
where x is the mole fraction of component 2 in 
the solid solution and 7 is a constant to be em- 
pirically determined. 

The special case where n=1 is known as 
Vegard’s law. This law has been found to hold 
true for a large number of solid solutions. Values 
of n=3 and n=8 have also been proposed. The 
latter value of m=8 has been proposed as a result 
of a theoretical study by Grimm and Herzfeld 
but has not been verified experimentally. The 
previously mentioned study of Havighurst, Mack, 
and Blake has shown that for values of 6 less 
than 6 percent solid solutions of the alkali 
halides at room temperature obey Vegard’s law 
rather closely. It is, however, hard to distinguish 
experimentally between the values n=3 and 
n=1 for these particular solid solutions. For 
simplicity in the theoretical derivation of the 
phase diagram, I have assumed the validity of 
Vegard’s law. 

By far the largest portion of the energy of an 
alkali halide is given by the Born evaluation of 
the electrostatic energy of the crystal lattice 
which is 

Uno = (NAe?2?/a9)(1—1/m),' 


where N = Avogadro’s law number, e= electronic 
charge, z=charge on ion, a9=distance from any 
ion to nearest neighbor, A = Madelung constant. 

The term in —1/m represents the repulsive 
energy of the lattice, where m is approximately 8. 
If we collect all the terms in the numerator of 
the above expression and call them /, we have 
Un =— M /ao. 

The Uoo represents by far the largest contri- 
bution to the energy of the ionic crystal lattice. 
Various other contributions such as van der 
Waal’s energy, polarization energy, etc., range 
from 1 to 5 percent as is seen from the values of 
Uy and U for the alkali halides listed in Table IT. 
U represents the total energy of the lattice® (not 
including thermal energy which is given by a 
Debye term). 

The heat of formation of these solid solutions 





* Uo is the electrostatic energy of the lattice at the 
absolute zero of temperature. 

®See Seitz, Modern Theory of Solids, Chapter 2. The 
more penetrating analysis of Mayer in which an exponential 
law for repulsion between the ions is used does not alter the 
results sufficiently to warrant its use in the approximations 
here used. 
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is in general endothermic, i.e., their formation 
absorbs heat. This is in accordance with the fact 
that solubility in the solid state increases with 
temperature. Also it is experimentally observed 
that the heat changes are greatest in the pairs 
KCl-NaCl, KBr-NaBr, KI-NaI whose mixed 
crystals are formed at high temperatures but are 
unstable at room temperature. Then comes KCI 
KI which has a large miscibility gap at room 
temperature, but at high temperature has a 
small miscibility gap. The heat of formation of 
the solid solution of KBr-KI is somewhat less; 
the miscibility which is complete just below 
melting shows at 25°C a gap between 34 and 93 
percent. Finally, the smallest heat of formation 
within the series is shown by KCI-KBr where 
solubility in the solid state is complete at room 
temperature.® 

If we assume the validity of Vegard’s law, we 
can make an approximate calculation of these 
heats of formation. The energy Uoo of an equi- 











TABLE II.{ 
—Uo —U calculated —U experimental 
(Sherman) (Mayer) (Sherman) 

Li F 240.1 242.0 

Li Cl 193.3 200.2 198.1 
Li Br 183.1 189.7 189.3 
Lil 170.7 176.2 181.1 
Na F 215.0 214.5 

Na Cl 180.4 184.4 182.8 
Na Br 171.7 176.5 173.3 
Na I 160.8 166.1 166.4 
K F 190.4 191.8 

K Cl 164.4 167.5 164.4 
K Br 157.8 161.7 156.2 
K I 149.0 153.1 151.5 
Rb F 181.8 184.2 

Rb Cl 158.9 162.7 160.5 
Rb Br 152.5 156.7 153.3 
Rb I 144.2 148.4 149.0 
Cs F 172.8 175.8 

Cs Cl 148.9 155.6 155.1 
Cs Br 143.5 150.0 148.6 
Cs I 136.1 142.7 145.3 











t See reference 5. 

Uo is calculated from the Born theory described above using values 
of m that are determined by use of experimental values of the com- 
pressibility and generally lie in the vicinity of 9. The lattice parameter 
used is that measured at room temperature and therefore U» is slightly 
different from Uoo because of thermal expansion. Mayer's more accurate 
calculations involve the use of an exponential term for the repulsive 
energy and also includes van der Waals energy. 


5 The foregoing remarks on heats of formation of these 
solid solutions are adapted from reference 1 as are the 
experimental data cited in Table III. 

















































190 








molal solid solution is given by 
Uoo= — M/([(ai+az)/2 ].7 


Let us take for example the reaction KCI+KBr 
= (solid solution KCI, KBr). (a; and az are here 
the lattice parameters of the pure 


2M MM 
aH=-———_+(—+—) 
(a; +a2)/2 a, do 


components which form the solid solution.) 











— 4a;d2.+42(a,+a2) +4;(a;+a2) 
aH =ar( ) 
(a, +42) (a1a2) 
M(a,—a2)*> M #& 


7 (ai +as)(aya2) 2a 10000 





approximately AH96? since 1//a~ 180,000 cal. 
(a is a kind of average between a; and az). 


We have therefore an approximate evaluation 
of the heats of formation of these solid solutions 
in terms of the deviation 


6 adga—aQ 
100 (as-+ay)/2 


In the above derivation we have used the fact 
that M/a=—Up> is on the average 180,000 
calories for the alkali halides. In Table III we 











TABLE III. 
AH (calculated) AH (observed) 
System 6 in calories in calories 

KCI-RbCI 4.53 166 203 

KBr-KCl 4.8 187 220 

NaCl-NaBr 5.5 236 aaa 
KBr-KI 7.1 389 (390) 
NalI-KI 9.1 646 (1250) 
NaBr-KBr 10.8 960 (1400) 
NaCl-KCl 11.5 1188 (2150) 











The numbers enclosed in parenthesis are unreliable because the 
equimolal solid solutions of the corresponding systems are unstable 
at room temperature. Since AH is obtained by measuring heats of 
solution in water this necessitates using quenched solid solutions which 
are of course thermodynamically unstable. 


7 Actually we should here use 


Uoo ao 1fao\” 9 6 
ae tals = (2) ]+gRe+RD(z) +0 
for the heat content of one mole of alkali halide crystal. 
But this is not very different from Uoo and we are using 
the latter as a good approximation. The heat of reaction 
is of course —AH. (9R0/8+RD(0/T) represents Debye 
energy.) @o represents the lattice parameter of absolute 
zero, a the lattice parameter at temperature in question. 
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have compiled the experimental and the calcu- 
lated values for the heats of formation of some 
of these solid solutions. It is only in the cases 
where Vegard’s law is obeyed and the miscibility 
is complete, i.e., cases where 6 is smaller than 
6 percent that we can expect heats of formation 
measured at room temperature to check the 
calculated values. 

In order to derive the phase curve for these 
solutions we consider the question: Suppose we 
mix in equimolal proportions two alkali halides, 
of lattice parameters a2 and a, with az2a,. If 
we allow the mixture to come to equilibrium the 
general result will be two phases, a solid solution 
of component 1 in component 2, and a solid 
solution of component 2 in component 1. It is a 
consequence of the assumptions and approxima- 
tions of this theory that the amount of com- 
ponent 2 in component 1 is equal to the amount 
of component 1 dissolved in component 2. This 
does not alter the general shape of the phase 
diagram and for simplicity we shall first derive 
the phase curve on this basis. We shall later 
discuss this question more fully from a theoretical 
standpoint and compare with experimental ob- 
servations. 

The total Gibbs free energy for one mole of 
alkali halide crystal is U+Aua+pV=F. U is the 
electrostatic energy of the crystal.* If we neglect 
expansion of the crystal, U= Uo. Aa represents 
the Debye expression for the Helmholtz free 
energy. PV is negligibly small below the melting 
point. 

Since a large portion of the free energy is com- 
prised by the term Upoo, and since the Debye 
temperatures and the extra electrostatic energies 
of expansion of the alkali halides are all of the 
same order of magnitude, we shall assume that 
insofar as it has any effect on the process of 
mixing F= Up. 

Now let us suppose that we mixed two alkali 
halide crystals (one mole of each). The free 
energy of the two solid solution phases formed is 


MM 
F=— —+—)+2RTIxIn x+(1—x)In(1—x)} 
a’ a 


where x is the mole fraction of component 1 in 
the one solid solution phase or the mole fraction 


8’ For a complete discussion see Slater, Introduction to 
Chemical Physics, Chaps. 23 and 13. 
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of component 2 in the other solid solution phase M a.—a, M az—ay 
and where a’ is the new lattice parameter of the a a 
. F ‘ moe P ay ay} do deo 
first solid solution phase and a” is the lattice sania 


parameter of the second solid solution phase. By de—a, 2 do—ay 2 
Vegard’s rule 1+(“—“)s in =| 


a, 


fon — 
a’ =a,+(d2—4a))x, +2RT7[In x—In (1—x)]=0. (2) 
a” =as+(a1—a2)x. Let us put 
Using these values ds— ay 
5.=———, 


M AM. ay 


Q:+(d2—a1)x d2+(a;—a2)x M 5; M be 
§a— — callie cape 


+2RT {x In x+(1—x) In (1—x)}. a; (1+6;x)* ae (1— dex)? 


In order to determine x we wish the free energy +2RT[In x—In (1—x)]. (3) 
to be a minimum. After carrying through a little... ~ 
titted S Since -6:26220.05; x<1 we can expand by 
algebra we can find x as a function of T and 6 ; e 
a , binomial theorem. 
the percent deviation of the two lattice param- 
eters. 


ete a 
a; 


0= M(a:—a)| —-— 


Qa;- do 


1 1 2M(a2—4a;)? 


M(a2—4a;) M(a,—a2) 





~ Fart (a2—a,)x 7 [a2+(ai—a2)x }? 2M (a2—a)* 


—_ ——x+2RT[In x—In (1—x)], (4) 
+2RT7(In x—In (1—x))=0. (1) ay3 : 
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M (ad2—4a)?(@2+a1) (=) 
= —— —}x 
1 1 


@179" a a 





2M dea 2 

-—| x+2RT7[In x—In (1—x)], (5) 
Qo Qo 

2M (d2—4a,)?/a*(1 — 2x) 


+2RT7[In x—In (1-—x)]=0. (6) 


(a is a kind of average between a; and a2.) For 
the alkali halides we have seen that the average 
value of Uoo is 180,000 cal. For the purpose of 
obtaining a general equation suitable for all 
alkali halides we set 


M/a~180,000 cal., R=2 cal./degree, 
6= 100(a2—a;)/a. 


Hence 98/7(1—2x)+1n x—In (1—x) =0. We let 


T 1-—x 
y=—, y=1=2+ /|in ( )} 
96° x 


om a 
VYmaximum > 2° 





We can make a table of y against x. See Table IV. 

The graph (Fig. 1) shows a plot of the mole 
fraction dissolved against 7/96. The charac- 
teristic feature of the curve is that it rises very 
steeply near the axis of zero concentration and 
then rounds off to reach a maximum beyond 
which miscibility is complete. The critical mixing 
temperature for various values of 6 is given in 
Table V. It is seen that this theory predicts 
complete miscibility at room temperature for 
values of 6 less than 8 percent as compared with 
the value of 5 or 6 percent that would be ex- 
pected on the basis of experiment. At higher 


TABLE IV. 


*0.01 0.02 0.05 0.10 0.20 0.30 040 0.50 
y 0.2175 0.2473 0.3004 0.3640 0.4329 0.4723 0.4938 0.5000 
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TABLE V. 
Critical mixing Critical mixing 
6 temperature 6 tem perature 
1 4.5° Abs. 10 450.0° Abs. 
2 18.0 11 544.5 
3 40.5 12 648.0 
4 72.0 13 760.5 
5 112.5 14 882.0 
6 162.0 15 1012.5 
7 220.5 16 1152.0 
8 288.0 17 1300.5 
9 364.5 | 18 1458.0 
| 19 1714.5 


| 
| 





temperatures the theory is in accord with the 
experimental data obtained so far. The initial 
steepness of the curve explains why many mix- 
tures completely miscible at high temperatures 
give at room temperature an x-ray photograph 
in which the two phases appear nearly pure. 

It is a consequence of the approximations of 
this theory that the mole fraction of component 
1 dissolved in component 2 is equal to the whole 
fraction of component 2 dissolved in component 
1. This can be seen by assuming that a certain 
amount x of component 1 is in solid solution 1 
and a certain amount y of component 2 is in 
solid solution 1. The composition of the other 
solid solution phase (solid solution 2) is then 
determined by the original amounts of the pure 
components. By equating the partial potentials 
of each component in the two phases we find 
that if we start out with equal amounts of the 
two pure components, x will be equal to y. The 
partial free energies are of course obtained by 
setting up the expressions for the total free 
energy of each phase according to the same 
approximations used above and differentiating 
with respect to the amount of component in the 
phase. The consequence that x=y is not borne 
out by experiment, but the general shape of the 
phase curve should not be radically wrong 
because of this error. 
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The Thermal Reaction Between Hydrogen and 
Oxygen at Higher Pressures 


O. OLDENBERG AND H. S. SOMMERS, JR. 


Research Laboratory of Physics, Harvard University, 
Cambridge Massachusetts 


January 27, 1942 


N recent experiments we found conditions under which 

the thermal hydrogen-oxygen reaction (p= 60 cm; T up 
to 560°C) is sufficiently reproducible so that the tempera- 
ture dependence of the rate could be measured. From the 
constant activation energy observed over a certain tem- 
perature range, we inferred that here the reaction proceeds 
in non-branched chains. In a recent letter to the editor, 
Dainton! concludes that the chains are branched. From an 
enlargement of our pressure-time curve representing a 
typical reaction, he concludes that at the beginning the 
first derivative of the rate (6?(Ap)/éf) is increasing and 
that the plot of log Ap against time is linear; applying 
Semenoff’s theory, he argues that in the first stages of the 
reaction branching of the chains is occurring and, in fact, 
more frequently than their termination. 

This conclusion seems questionable for the following 
reasons. The beginning of our curve analyzed by Dainton 
is not reproducible enough to apply Semenoff’s criterion; 
in many cases 6?(Ap)/éf? seems to begin with its largest 
value; in still more cases (old vessel) the reaction begins 
with nearly the maximum rate. The reproducibility of the 
maximum rate (or rather log rate), however, is evident from 
our curves (reference 1, Figs. 2-5). 

Moreover, immediately after high evacuation near the 
annealing temperature of Pyrex and subsequent admission 
of almost one atmosphere of gases, one should expect 
rapid changes of the surface by adsorption superimposed 
on the effects considered by Semenoff. This difficulty may 
be less pronounced in the low pressure explosion (Kowalsky, 
pressure a few mm, reaction lasting a small fraction of a 
second) ; this reaction, which is quantitatively described by 
Semenoff’s theory, takes place under conditions so widely 
different from those at high pressure that Hinshelwood and 
Williamson even postulated a mechanism of a different 
nature.* Whether the high pressure reaction should actually 
be interpreted by Semenoff’s theory can possibly be deter- 
mined by a photochemical investigation; if the induction 
period is correlated with the time required for the develop- 
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ment of a chain, a corresponding decay period should be 
observed after the interruption of the illumination. Taylor 
and Salley, in their photochemical investigation, mention 
the induction period but not such a decay period. 

There is another difficulty in Dainton’s interpretation. 
He assumes that in the beginning of the reaction the 
branching exceeds the breaking and the explosion is pre- 
vented by the change of concentrations occurring during 
the development of the first chains. He argues that the 
maximum rate occurs when @¢=0(6—8=0) and that in 
this region the chains are of constant finite length. While 
we are unable to find a proof for this general theorem, it 
seems to us that the condition for the maximum rate would 
largely depend on the detail of the reaction. In the case 
thoroughly treated by Semenoff,* “the velocity of the 
reaction begins to decrease long before 6—8 becomes equal 
to 0.” 

We applied the conventional procedure (Hinshelwood 
and Thompson, Taylor and Salley, and Prettre), measuring 
the nearly constant rate after the induction period. This 
seems justified in the limiting case of an initial rate which 
nearly equals and rapidly reaches the maximum rate (many 
observations in Pyrex, Prettre’s observations in the KCl 
covered vessel); here one should be able to derive an 
activation energy from the temperature coefficient of the 
maximum rate. The same procedure leads to well-defined 
values of the activation energy for various kinds of surface 
yielding short or long induction periods. Therefore, we 
consider our interpretation of our curves most plausible, 
although not inevitable. This procedure does not involve 
an assumption regarding the nature of the induction period. 
In order to explore the effects of a pretreatment of the walls 
and of impurities on this period, we started experiments a 
year ago which are now interrupted. 

1F. S. Dainton, J. Chem. Phys. 9, 826 (1941); see F. S. Dainton and 
R. G. W. Norrish, Proc. Roy. Soc. A177, 393 (1941). For further 
references see Dainton’s paper and our paper, J. Chem. Phys. 9, 432 
CN. Hinshelwood and A. T. Williamson, The Reaction between 


im and Oxygen (Clarendon Press, Oxford, 1934), p. 52. 
. Semenoff, Chemical Kinetics and Chain Reactions, p. 62. 





The Infra-Red Absorption of CO" at 
4.66 Microns 


R. T. LAGEMANN* 


Mendenhall Laboratory of Physics, The Ohio State University, 
Columbus, Ohio 


January 23, 1942 


FTER taking the necessary experimental precautions, 

it has been found possible to find rotational lines due 

to the isotopic carbon monoxide molecule C“O'. These 

lines are found interspersed among the lines of the funda- 
mental band of CO" at 4.66 microns. 

A portion of the absorption curve is shown in Fig. 1. 
The more intense lines are due to the C"O"* molecule. 
The weaker lines are associated with the C“O" molecule, 
and comprise part of the R branch of that molecule’s 
fundamental band. At lower frequencies than are shown 
here the isotopic lines become relatively more intense. As 
is to be expected, the line spacing for the heavier isotope 
is less than that for the lighter one. 
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The grating used was an echelette replica of 2887 lines 
per inch, ruled for the region 9 to 12 mu and employed in 
the second order. Readings of percent absorption were 
made at intervals of 5 seconds of arc, equivalent to about 
0.13 cm. The slit width was equivalent to 0.22 cm™. 

Measurements are soon to be made on the overtone band 
and a complete report will be made at a later date. 


* Emory University, Emory University, Georgia. 





The Vibrational Structure of Electronic Transi- 
tions for Some Complex Ions 


MELVIN L. SCHULTZ 
Department of Chemistry, Illinois Institute of Technology, Chicago, Illinois 
February 9, 1942 


OR several compounds of some of the transition ele 
ments, absorption spectra consisting of groups of 
approximately equidistant bands have been reported. The 
positions of the bands and the intervals between them 
depend only upon the metal ion and upon the identity of 
the complex of which the metal ion is a part. These spectra 
may be interpreted as originating in electronic transitions 
in which the vibrational structure characteristic of the 
complex ion present is resolved. The experimental basis 
for this interpretation is summarized in the following 
paragraphs. 

Dreisch and Trommer! have found in dilute aqueous 
solutions of CoClz:, CoBr2 and Col, an absorption band at 
1.23 consisting of six maxima with an average separation 
of 248 cm. This band may be attributed to the Co(OH2).6** 
ion. In aqueous solutions of NiCl2, NiSO, and Ni(NOs)> 
they have found at 1.17% a band characteristic of the 
Ni(OH2).6** ion. This band is resolved into four maxima 
with a separation of 264 cm™. In ammoniacal solutions of 
Ni(NHs)6Cl2 the band is shifted to 0.940u and the separa- 
tion becomes 390 cm, the number of maxima remaining 
the same. For solutions of a number of ferrous salts, 
Dreisch and Kallscheuer? have reported a series of twelve or 
thirteen bands extending from about 0.74 to about 1.27y. 
The intervals separating the five bands of highest frequency 
vary considerably from salt to salt. The separations of the 
next seven bands are approximately the same for all of the 
salts measured. This constant separation, 324 cm™, is then 
characteristic of the Fe(OH2)¢** ion. 

The absorption spectra, measured at liquid air tempera- 
tures, of solid hydrates of some of the above salts have been 





reported.’ These spectra consist of many bands, from about 
20 cm™ to about 150 cm™ in width, extending throughout 
the visible region. Analysis indicates that these spectra 
may, in part, be interpreted in the same manner as those 
from the solutions. With CoCl,-6H.O, the average wave 
number difference for twelve pairs of bands is 226 cm“. 
With CoBr:.-6H.0, the difference for ten pairs of bands is 
220 cm™. For NiCl.-6H,0, the difference for nine pairs of 
bands is 237 cm. For Co(OH2)6**, in passing from solu- 
tion to the crystal, there is a decrease of 25 cm™ in the 
separation of the vibrational bands. For Ni(OH2).**, there 
is a corresponding decrease of 27 cm™. 

The interpretation presented here is consistent with the 
information obtained from the Raman effect of amino 
complex ions.4 Cu(NH3)4**, Zn(NH3)6**, and Cd(NHs3).6*" 
give Raman displacements of the same order of magnitude 
as the separations of the vibrational levels of Ni(NHs3)«6** 
and the aquo complexes (see Table I). If one makes the 


TABLE I. 














Force 

Vibration frequency constant 

Complex ion cm~! dynes/cm 

Cu(NHs3)4** (solution) 415 (Raman) 1.16 X10" 
Cd(NHs3)6** (solution) 340 (Raman) 0.85 
Zn(NHa)6** (solution) 423 (Raman) 1.22 
Ni(NHs3)«6** (solution) 390 (Absorption) 1.01 
Ni(OH2)6** (solution) 264 (Absorption) 0.52 
Ni(OH2)6** (crystal) 237 (Absorption) 0.42 
Co(OH2)s** (solution) 248 (Absorption) 0.46 
Co(OH:2)«** (crystal) 223 (Absorption) 0.37 
Fe(OH:2)«** (solution) 324 (Absorption) 0.77 
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approximation that the vibration of the complex ion may 
be represented by a simple harmonic oscillator, the force 
constant of the bond between the central metal ion and the 
coordinated group may be estimated from the vibration 
frequency. The values so obtained are listed in the table. 

Some evidence for anharmonicity is obtained from a 
series of five bands for solid CoCl.-6H20.4 The wave 
numbers of these bands may be represented with an 
accuracy of one wave number by: 


v=24504+ 222n—1.6n? 


From the constants of this equation, the heat of dissociation 
of Co(OH:2)¢** in crystals of CoCl,-6H:0 is estimated to be 
22,000 cal. The heat of the reaction: 


CoX2:6H20(s)+CoX 2(s) +6H20(/) 


determined from thermochemical data (I.C.T.) is 21,300 
cal. for X =Cl, 19,800 cal. for X =Br, and 16,700 cal. for 
X = NOs. The decrease in the thermochemical value for the 
bromide as compared with the chloride parallels the de- 
crease in the average separations of the vibrational levels, 
226 cm= for the chloride as compared with 220 cm™ 
for the bromide. The nitrate has not yet been studied 
spectroscopically. 


(n=0, ---, 4). 


Th. Dreisch and W. Trommer, Zeits. f. physik. Chemie B37, 37 
(1937). ; 

2 Th. Dreisch and O. Kallscheuer, Zeits. f. physik. Chemie B45, 19 
(1939). 

3 Gieleszen, Ann. d. Physik [5] 22, 537 (1935). 


41. Damaschun, Zeits. f. physik. Chemie B16, 81 (1932). 








